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1 Introduction 

Field  investigation  in  Northern  Mississippi  has  demonstrated  that  the  impact  of  LWD 
Jams  changes  downstream  through  the  channel  network.  This  downstream  change  is 
conceptualized  in  the  Debris  Jam  Classification  Model  (see  Wallerstein  et  ah,  1997).  It  is 
very  difficult,  however,  to  precisely  determine  which  components  of  channel 
geomorphology  are  attributable  to  LWD  influences,  and  which  are  caused  by  other,  more 
dominant  influences  such  as  channel  incision  processes.  It  is  also  difflcult  to  determine 
cause  and  effect  when  assessing  the  true  geomorphic  impact  of  LWD  in  streams. 

In  order  to  overcome  the  problem  of  interpretation,  a flume  model  offers  an  excellent  tool 
for  more  accurately  determining  the  hydraulic  forces  exerted  upon  LWD  elements  and 
the  impact  of  energy  dissipation  around  the  elements  upon  the  channel  boundary.  In 
order  to  obtain  meaningful  experimental  results  it  was  necessary  to  create  a scale  model 
of  a field  site  where  survey  of  LWD-induced  channel  change  had  been  conducted,  so  that 
model  results  could  be  assessed  against  field  data.  This  report  presents  an  overview  of 
research  carried  out  by  others  on  the  hydraulic  effects  of  LWD,  the  theory  of  scale 
physical  modeling,  the  derivation  of  parameters  for  the  model  created,  a discussion  of  the 
parameters  measured  and  their  relative  importance,  and  the  presentation  and  analysis  of 
the  results  obtained. 
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2 Literature  Review  of  Hydraulic  Impacts  of  LWD 

Previous  studies  have  investigated  the  effect  of  LWD  on  runoff  hydrographs,  velocity 
distributions,  water  surface  profile,  and  channel  roughness. 

2.1  Effect  of  LWD  on  stage/discharge  relationships,  the  hydrograph  and  flood 
frequency 

LWD  clearly  influences  the  direction  and  magnitude  of  flow  currents  within  stream  flow,  but 
few  data  have  been  documented  in  the  literature.  Swanson  and  Leinkaemper  (1978) 
produced  detailed  maps  of  debris  jams  indicating  flow  direction.  Smith  and  Shields  (1992) 
reported  that  the  removal  of  LWD  from  a river  18  to  23  m wide  and  3.5  to  4.5  m deep 
produced  more  uniform  flow,  with  less  of  the  channel  occupied  by  eddies  or  regions  of 
reduced  velocity. 

LWD  is  often  removed  from  channels  because  it  is  believed  this  will  significantly  reduce 
channel  roughness  which  will,  in  turn,  produce  a higher  mean  velocity  and  increase  in-bank 
channel  flow  capacity.  There  is  some  evidence  to  support  this  assumption.  For  example. 
Smith  and  Shields  (1992)  measured  the  mean  velocity  in  two  cleared  reaches  of  the  Obion 
River,  Tennessee  to  be  0.4  m/s  and  0.34  m/s  respectively,  while  in  an  uncleared  reach  of  the 
same  river  the  mean  velocity  was  0.27  m/s.  MacDonald  and  Keller  (1987)  found  that  there 
was  a local  increase  in  velocity  by  up  to  250%  as  a result  of  LWD  removal  and  decrease  in 
sinuosity  of  the  low-flow  thalweg.  According  to  Gippel  et  al.  (1992),  the  Murray-Darling 
Basin  Commission  calculated  a theoretical  reduction  in  water  level  of  0.3  to  0.4  m after  the 
removal  of  approximately  200  debris  jams  per  kilometer.  However,  later  analysis  of  flow 
records  indicated  a reduction  of  only  0.2  in.  There  should  be  a statistical  reduction  in  the 
magnitude  and  frequency  of  overbank  flooding  where  debris  is  removed  from  a channel 
because  of  the  increased  channel  capacity.  Bodron  (1994)  used  a dynamic  routing  model  to 
demonstrate  changes  in  both  stage  and  duration  of  flood  events  before  and  after  LWD 
removal,  using  values  of  Manning's  n calculated  in  the  study  by  Smith  and  Shields  (1992)  at 
South  Fork  Obion  River,  West  Tennessee.  Despite  the  fact  that  the  increase  in  channel 
cross-sectional  area  due  to  LWD  removal  was  ignored,  small  reductions  in  flood  height  and 
duration  were  calculated  based  solely  on  the  change  in  Manning's  n.  Bodron  (1994)  also 
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noted  that  flood  stage  would  be  reduced  further  if  sediment  accumulations  at  each  Jam  site 
had  been  removed.  However,  according  to  Gippel  et  al.  ( 1992),  many  claims  regarding  flood 
routing  lack  supporting  data.  It  is  even  possible  that  LWD  removal  might  increase 
downstream  flood  peaks  because  in  the  smoother  channel  the  flood  wave  is  less  attenuated. 
Gregory  et  al.  (1985)  found  that  LWD  ponds  water  which  results  in  an  increase  in  water 
depth  and  a decrease  in  velocity  which,  at  low  flows,  influences  travel  time  significantly.  At 
high  flows,  however,  the  ponding  effect  of  LWD  is  drowned  out.  Shields  and  Nunnally 
( 1984)  noted  that  because  large  accumulations  of  LWD  have  a damming  effect  on  the  flow 
which  locally  elevates  the  base  level  they  can  be  treated  as  geometric  elements  within  the 
channel,  rather  than  simply  as  roughness  elements,  in  backwater  profile  computations. 

2.2  Modeling  the  impact  of  LWD  on  flow  afflux 

Most  studies  of  resistance  to  flow  in  rivers  have  concentrated  on  small-scale  roughness,  such 
as  grain  or  bedform  roughness  where  the  size  of  the  roughness  element  is  small  compared  to 
the  flow  depth.  LWD,  on  the  other  hand,  represents  a form  of  large-scale  roughness  (Petryk 
and  Bosmajian,  1975).  Flow  conditions  associated  with  the  presence  of  LWD  in  streams 
varies  from  sub-critical  to  super-critical  depending  on  the  dimensions  of  the  LWD  and  the 
depth  of  water. 

Gippel  et  al.  (1992)  used  the  momentum  equation  to  detennine  the  hydraulic  effect  of  LWD 
in  terms  of  an  afflux  or  backwater  effect.  If  flow  is  subcritical  (Froude  number  < 1 ),  then 
apart  from  loeal  disturbance  of  the  veloeity  profile,  LWD  only  has  an  influence  in  the 
upstream  direetion.  Quantifying  the  backwater  effects  is  problematical  because  of  the 
practical  difficulties  of  directly  measuring  the  afflux  at  debris  jams.  An  alternative  to  direet 
measurement  is  prediction  on  the  basis  of  a known  relationship  between  afflux  and  more 
easily  measured  parameters.  Gippel  et  al.  (1992)  used  the  results  of  a laboratory  hydraulic 
study  to  develop  a method  of  determining  the  afflux  caused  by  LWD  (see  Figure  1 ). 

Gippel  et  al.  (1992)  propose  the  use  of  the  following  equation  to  calculate  afflux: 
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(Fr'  -l)  + V(Fr'  -l)  +3CoBFr' 


Ah  = 


where.  Ah  = afflux  = h|  - h2  (m),  and  the  drag  coefficient  (Cd)  is  given  by 


C D 


Fd 

K pv“A 


(2) 


where,  Fd  = drag  force  (N),  p = density  of  water  (1000  kg/m^),  v = mean  approach  velocity 
to  object  (m/s),  and  Ad  = area  of  the  LWD  element  projected  in  the  How  direction  (m“). 
The  Froude  number  (Fd  is  defined  as: 


Fr  = ^X^  (3) 

where  V2  = mean  velocity  at  section  downstream  of  the  LWD  (m/s),  and  h2  = water  depth 
downstream  of  LWD  (m) 


The  blockage  ratio  (B)  is  defined  as: 

B = Ad/A  (4) 


where  A = w (channel  width)  x hi  = cross  sectional  area  of  flow  (m"). 


Figure  1:  Definition  sketch  of  LWD  model  by  Gippel  et  al.  (1992) 
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Thus,  the  afflux  depends  on  Fr,  Cd  and  B.  The  Froude  number  can  be  calculated  from  direct 
measurement  or  from  flow  records.  B can  be  found  from  a field  survey. 

The  remaining  problem  centers  on  selecting  an  appropriate  drag  coefficient.  The  drag 
characteristics  of  a cylinder  in  infinite  flow  are  well  known  (Petryk  and  Bosmajian,  1975). 
Less  is  known  about  drag  on  cylinders  within  boundaries  where  the  “blockage  effecf’  is 
significant  and  the  drag  coefficient  is  consequently  increased.  Gippel  et  al  ( 1992)  conducted 
experiments  on  LWD  models  to  determine  drag  force,  using  a towing  carriage  and  water 
tunnel.  Froude  number,  LWD  length  to  diameter  ratio,  and  LWD  distance  from  the  bed  all 
affected  drag  eoefflcient.  but  were  much  less  important  than  the  blockage  effect,  angle  of 
orientation  to  the  flow,  and  the  shielding  effect  (i.e.,  one  piece  of  LWD  behind  another).  A 
suitable  drag  coefficient  (C'n)  for  the  LWD  in  question  can  be  selected  from  their 
experimental  results  (Gippel  et  al.  1992,  figures  3.8  or  3.12)  on  the  basis  of  its  overall  shape 
and  angle  of  orientation.  The  drag  coefficient  should  then  be  adjusted  for  the  blockage 
effect,  which  can  be  calculated  using  the  following  equation  developed  by  Gippel  et  al. 
(1992)  using  their  empirical  data  from  flume  studies; 

Co  = Cd(1-B)-^  (5) 

where  Cd  - adjusted  drag  coefficient,  and  C'd  = drag  coefficient  in  infinite  flow.  These 
data  are  then  substituted  into  equation  I to  calculate  the  afflux. 

Predicted  and  measured  afflux  values  resulting  from  the  flume  study  were  very  closely 
correlated,  and  Gippel  et  al.  (1992)  concluded  that  the  flume  conditions  did  not  seriously 
violate  any  of  the  assumptions  in  equation  1 . The  proposed  method  of  afflux  estimation  was 
then  applied  to  data  collected  from  the  Thomson  River,  Vietoria  and  revealed  that  de- 
snagging  would  produce  a reduction  in  stage  of  only  0.01  m at  bankfull  flow. 

In  conclusion  the  authors  believe  this  method  of  backwater  or  afflux  calculation  due  to 
individual  LWD  fonnations  could  be  used  as  a tool  to  help  determine  whether  the  afflux 
reduction  due  to  LWD  removal  would  have  a significant,  positive  impact.  Conversely,  this 
method  could  determine  whether  LWD  could  be  left  in  place  and  perhaps,  re-orientated. 
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lopped  or  even  re-introduced  where  sympathetic  rehabilitation  management  is  desirable, 
without  significant  effect  on  high  in-bank  stages. 


2.3  Effect  of  LWD  on  channel  roughness 

The  Manning  equation  is  based  on  a resistance  coefficient  that  represents  all  sources  of 
roughness  in  the  channel  and  is  defined  as: 

2 I 

R 'S  2 

n = (6) 

V 

where  n = Manning's  roughness  coefficient.  R = hydraulic  radius  (m),  and  S = energy  slope. 
This  equation  is  widely  used  by  river  engineers  who  select  values  of  n from  personal 
experience,  tables  in  Chow  (1959)  or  photographs  in  Barnes  (1967).  The  range  of  n values 
in  stream  channels  without  LWD  is  between  0.025  and  0.15.  For  streams  heavily  conjested 
with  LWD,  less  than  30  m wide,  n ranges  from  0.075  to  0.15.  Irregular  and  rough  reaches  of 
large  streams  have  values  of  n from  0.035  to  0. 1 0. 


Manning’s  n has  an  empirical  derivation,  however,  while  the  Darcy-Weisbach  flow 
resistance  equation  has  the  advantage  of  being  theoretically  based  (Richards,  1982).  The 
Darcy-Weisbach  f is  defined  as: 

8gRS 


f 


(7) 


where  f=  Darcy-Weisbach  friction  factor  and  g = acceleration  due  to  gravity  (9.81  m/s“). 


The  effect  of  LWD  on  flow  resistance  varies  as  a function  of  relative  flow  depth.  Bevan  et  al. 
(1979)  found  that  when  LWD  is  large  in  relation  to  flow  depth  the  roughness  coefficient  is 
extremely  high  (Manning's  n>l).  As  LWD  is  submerged  it  e.xerts  less  influence  on  flow 
resistance.  Smith  and  Shields  (1992)  measured  a large  decrease  in  f as  discharge  increased. 
They  also  observed  that  friction  factors  for  cleared  and  uncleared  reaches  converged  at  high 
in-bank  flows.  Indirect  evidence  to  support  these  findings  is  provided  by  investigations  of 
downstream  hydraulic  geometry  which  show  that  roughness  generally  decreases  as  channel 
size  increases  (Wolman,  1955).  Petryk  and  Bosmajian  (1975)  derived  the  following  equation 
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to  predict  n as  a function  of  density  of  vegetation  in  the  channel,  hydraulic  radius,  n due  to 
boundary  roughness  and  a vegetation  drag  coefficient  (in  US  units): 


where  nb  = Manning's  boundary  roughness  coefficient  excluding  the  effect  of  vegetation,  Cd 
= drag  coefficient  for  vegetation.  Ad,  = projected  area  of  the  /th  plant  in  the  streaniwise 


perimeter  (m).  In  this  fomiula  the  expression  CdSAd,/AL  represents  the  density  of 
vegetation  in  the  channel.  A problem  with  this  formula  lies  in  selecting  a value  for  the  drag 
coefficient,  Cd-  Petryk  and  Bosmajian  (1975)  assumed  a value  of  I,  but  as  mentioned 
previously  this  strictly  only  applies  to  cylinders  in  infinite  fiow.  In  streams,  interference 
from  nearby  obstructions  and  the  effect  of  blockage  on  the  drag  coefficient  must  also  be 
considered.  These  effects  are  discussed  further  in  section  6.3. 

Manning's  equation  is  also  inapplicable  in  situations  where  there  is  a high  degree  of 
obstruction  in  the  channel,  particularly  where  n>l.  This  is  because  n was  empirically 
derived  to  describe  open  channel  situations  with  fully  turbulent  flow  where  friction  is 
controlled  primarily  by  skin  friction  at  the  channel  boundary.  The  equation  attaches 
significance  to  the  hydraulic  radius  which  may  be  irrelevant  if  the  channel  is  heavily  choked 
with  LWD.  It  was  therefore  considered  inappropriate  to  use  Petryk  and  Bosmajian's  (1975) 
equation  to  determine  average  friction  factor  values  for  debris  elements  in  the  study  reaches. 

Smith  and  Shields  (1992)  studied  the  effects  of  varying  levels  of  LWD  density  on  the 
physical  aquatic  habitat  of  South  Fork  Obion  River,  Tennessee,  USA.  Two  secondary 
objectives  in  their  study  were  to  develop  and  demonstrate  a method  for  quantifying  LWD  in 
a given  reach  and  to  relate  the  quantity  of  LWD  to  reach  hydraulics.  They  used  a theoretical 
approach  similar  to  that  described  by  Petryk  and  Bosmajian  (1975)  to  calculate  the  effect  of 
LWD  on  channel  roughness,  but  based  their  analysis  around  deriving  f rather  than  n.  The 
LWD  density  in  a reach  was  calculated  using  the  following  fomuila.  (Smith  and  Shields, 
1992): 


(8) 


direction  (m^),  L = length  of  the  channel  reach  being  considered  (m),  and  P = channel  wetted 
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where  DA  = LWD  density  (m‘),  n = total  number  of  LWD  formations  in  the  reach,  Fb,  = 
formation  type  weighting  factor  for  /th  formation  type,  N,  k = number  of  type  j LWD 
formations  in  klh  width  category,  and  Fwk  = weighting  factor  based  on  LWD  formation  width 
category.  Factors  Fb,  N,  Fw,  k and  j are  defined  on  the  LWD  survey  form  devised  by  Smith 
and  Shields  which  is  presented  in  the  Appendix. 


In  a channel  reach  where  LWD  plays  a major  role  in  flow  resistance,  total  resistance  can  be 
expressed  as  (Smith  and  Shields,  1992): 

ft=tb+Lebns  (10) 

where  tj  = total  Darcy-Weisbach  friction  factor,  fb  = boundary  friction  factor  excluding  LWD 
effects,  and  fdedns  ^ friction  factor  due  to  LWD.  They  then  define  total  head  loss  over  a reach 
as  the  sum  of  a boundary  friction  loss  and  a LWD  blockage  loss  as  follows; 


2g 


(11) 


where  h|=  total  head  loss  (m),  and  Kd  = dimensionless  loss  coefficient  (dependent  upon 
LWD  density).  The  energy  gradient  (S)  is  calculated  using  a total  friction  factor  from  the 
Darcy-Weisbach  equation  (equation  7). 


By  substituting  equation  7 for  S into  equation  1 1,  Smith  and  Shields  (1992)  obtained: 


f , = f b + 


4RKd 

L 


(12) 


thus: 


4RKd 

fdebns  — (13) 

The  ratio  Kd/L  is  then  expressed  in  terms  of  the  LWD  density  as  (Smith  and  Shields,  1992): 
Kd/L  = DA  (14) 


Smith  and  Shields  (1992)  calculated  values  for  fb  using  curves  developed  by  Alain  and 
Kennedy  (1969),  hydraulic  parameters  detennined  from  dye  tracer  tests  in  the  LWD  reaches 
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(which  provide  direct  discharge  and  velocity  estimates  (Richards  1982)),  and  the  median  bed 
grain  size  determined  from  sieve  analysis.  Values  for  fdebns  were  then  calculated  using 
equations  8,  14  and  15.  They  then  compared  computed  values  of  ft  with  values  measured 
using  dye  tests. 

The  results  of  Smith  and  Shields  (1992)  study  showed  a positive  correlation  between  the 
measured  and  computed  friction  factors.  However,  they  reeognized  that  considerable 
refinement  and  site-specific  adaptation  may  be  required,  and  that  the  method  does  not 
account  for  local  energy  loss  because  of  bends  or  flow  expansion  and  contraction  at  bridges, 
debris  dams,  or  riffles.  I’he  method  does  have  a sound  theoretical  basis,  however,  and  could 
be  usefully  employed  in  future  research  into  the  hydraulics  effect  of  LWD.  This  analysis 
was  further  refined  by  Shields  and  Gippel  ( 1995)  who  produced  a table  of  fdebns  values  for 
LWD  jams  on  the  Obion  River  (Tennessee)  and  the  Tumut  River  (Australia).  Their  results 
are  used  to  compare  and  validate  fdebns  values  obtained  from  the  flume  model  (see  below). 

2.4  Flume  techniques 

Only  two  additional  papers  have  been  found  concerning  experimental  work  on  LWD 
dynamies  from  a eomprehensive  literature  survey.  Young  (1991)  conducted  flume 
studies  to  investigate  the  impact  of  debris  obstruction  on  flow  afflux,  while  Cherry  and 
Beschta  (1986)  used  a flume  with  an  erodible  bed  to  determine  the  effect  of  debris 
orientation  and  height  above  the  bed  on  the  depth  and  extent  of  bed  scour.  Neither  of 
these  studies  made  any  attempt  to  scale  flume  parameters  from  a field  prototype,  so 
results  obtained  in  these  studies  are  useful  only  to  deseribe  the  general  nature  of  the 
processes  involved. 
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3 Similarity  in  Model  Scaling 


3.1  Basic  principles 

Model  scaling  depends  upon  dimensional  analysis  whereby  the  number  and  complexity 
of  experimental  variables  which  affect  a given  physical  phenomenon  are  reduced  or 
compacted  (White,  1994).  There  are  only  four  basic  dimensions  which  need  be 
considered  in  any  mechanics  problem.  These  are  Mass  (M),  Tength  (L),  Time  (T).  and 
Temperature  (9).  This  is  the  MLT0  system  (note  that  Force  (F)  sometimes  replaces  mass, 
to  give  the  FLT0  system).  Dimensional  analysis  makes  use  of  the  MLT0  system  to  create 
scaling  laws  which  enable  one  to  construct  models  which  faithfully  replicate  the  relative 
magnitudes  of  the  most  important  forces  that  are  present  in  the  Teal  world'  prototype. 

Before  scaling  a model  one  must  consider  which  parameters  are  important  in  that  system, 
i.e.  which  variables  one  wishes  to  measure.  This  choice  introduces  a degree  of  ambiguity 
into  the  method,  so  the  engineer  must  think  carefully  before  selecting  a group  of 
variables. 

In  the  case  of  an  open  channel  flow  tluid  mechanics  problem  one  is  essentially  interested 
in  pressure  change  (Ap)  through  a fluid  over  space  and  time  (Henderson,  1966).  If  flow 
is  sub-sonic  pressure  is  dependent  upon: 

Ap  = f (v,  L,  p,  )u,  c,  g)  (15) 

where  L = a length  scale,  usually  flow  depth  (d)  (m),  p = fluid  density  (approx.  1000 
kg/m^  at  20°C;  Richards,  1982),  p = tluid  viscosity  (approx.  0.001  m“/s  at  20°C; 
Richards,  1982),  and  a = fluid  surface  tension  (approx.  0.073  N/m  for  a water-air 
interface;  Chadwick  and  Morfett,  1991) 

At  the  core  of  dimensional  analysis  is  the  Buckingham  pi  (IT)  theorem  (see  White,  1994; 
Roberson  and  Crowe,  1993).  This  theorem  states  that  for  any  group  of  variables  selected 
there  are  n-k  dimensionless  groups  that  govern  the  system,  where  n is  the  number  of 
variables  and  k is  the  number  of  dimensions  of  importance.  If  we  assume  that 
temperature  is  not  of  great  importance  in  the  case  of  a free-surface  fluid  flow  problem,  k 
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= 3 (MLT).  In  most  modeling  studies  the  seven  variables  listed  in  equation  15  are  usually 
considered  to  be  of  the  greatest  importance.  This  suggests  that  there  are  7 - 3 = 4 
important  dimensionless  groups  in  the  system.  The  variables  listed  in  equation  15  have 
the  following  MLT  dimensions: 


Quantity 

Symbol 

MLT  UNITS 

Pressure 

P 

M/LT- 

Velocity 

V 

ITT 

Length 

L 

L 

Density 

P 

M/0 

Viscosity 

P 

M/LT 

Surface  Tension 

a 

M/T- 

Gravity 

a 

L/T- 

Table  1:  Dimensions  for  important  fluid  mechanics  variables 


These  variables  are  then  compacted  into  dimensionless  groups  according  to  their 
dimensional  properties.  The  groups  are  not  derived  here  as  the  method  is  discussed 
thoroughly  elsewhere  (see  Henderson,  1966,  Yalin,  1971,  White,  1994).  This  technique 
produces  the  following  four  dimensionless  groups  (Roberson  and  Crowe,  1993): 


'2  pW 


= / 


pvL  p\' 


\ ^ 


a 


or  C = /(Re,  We,Fr) 


(16) 


where  Cp  = Pressure  coefficient.  Re  = Reynolds  number,  Fr  = Fronde  Number,  and  We  = 
Weber  number.  Because  pressure  distribution  in  a fluid  is  an  unknown  function  of  Re, 
Fr,  and  We,  model  studies  can  only  be  interpreted  accurately  if  Re,  Fr,  and  We  are 
severally  given  the  same  value  in  the  model  as  in  the  prototype. 


If  a model  is  to  have  complete  similarity  with  a prototype  there  must  be  geometric 
similarity  (the  model  looks  the  same  as  the  prototype),  dynamic  similarity  (model  forces 
are  scaled  proportionally  to  the  prototype),  and  kinematic  similarity  (where  flow  patterns 
in  the  model  map  onto  those  of  the  prototype).  This  requires  that  all  velocities  and 
accelerations  must  have  the  same  prototype  to  model  ratio,  such  that: 

(Va)p  _ (Vb)p 
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where  subscripts  a = at  section  'a',  b = at  section  ‘b\  p = prototype,  and  m = model. 


3.2  Secondary  scale  ratios 

Building  a model  requires  that  secondary  scale  ratios  be  developed  to  translate  prototype 
quantities,  such  as  discharge  and  slope,  into  corresponding  model  values.  The  influence 
of  the  Froude  number  is  dominant  if  the  channel  has  a free  surface,  so  it  is  of  critical 
importance  that  Fr  has  the  same  value  in  the  model  as  in  the  prototype.  Thus  (Roberson 
and  Crowe,  1993): 


Having  decided  that  the  ratio  of  the  Froude  number  between  model  and  prototy  pe  must 
equal  unity,  we  can  derive  a set  of  scale  ratios  for  MLT  as  follows. 


First,  given  that  g is  a universal  constant,  regardless  of  scale,  the  Froude  number  equation 
can  be  rearranged  to  give  a length  scale  ratio  (Lr)  in  terms  of  flow  velocity: 


V2 


L. 


or 


Lr  = Vr" 


(19) 


where  the  subscript  r refers  to  the  ratio  of  prototype  to  model.  Second,  the  time  scale 
ratio  (Tr),  in  order  to  remain  faithful  to  Froude  scaling,  must  also  be  deduced  as  a 
function  of  the  velocity  scale  and  length  scale.  Velocity  has  MLT  units  of  L/T,  and  the 
reciprocal  is  therefore  that  T = Lv"',  thus  (Yalin,  1971): 

Tr  = LrVr-'=Lr"’  (20) 


Finally,  mass  is  defined  by  the  density  of  an  object  or  fluid,  multiplied  by  its  volume,  so 
mass  scale  ratio  (Mr),  is  equal  to  (Yalin,  1971 ): 

Mr^PrLr^  (21) 


We  have  now  derived  a set  of  scale  ratios  for  MLT  based  upon  unity  of  the  ratio  of 
Froude  number  between  model  and  prototype,  so  all  other  important  scale  ratios  can  now 
be  derived  similarly  using  equations  19,  20  and  21.  The  important  scale  ratios  are  shown 
in  Table  2. 


.)  V'  * . .1 


f. 


• ■ I 


!> 


I 


'.1 


( 


1\  f 

'tl  ^ i') 


I - Ji ' 


/'  I 


ii 


i ^ 

g thji 


Variable  Ratio 

Symbol 

Scale  Ratio 

Mass 

Mr 

prLr^ 

Length 

Lr 

Lr 

Time 

Tr 

I “ 1 

LrVr  - 

Lr"^ 

Velocity 

Vr 

Lr"^ 

Discharge 

Qr 

VrLf  “ 

Lr=^ 

Force 

Fr 

MrLrTr 

T 3 

PrLf 

Pressure 

Pr 

II 

PrLr 

Table  2:  Important  scale  ratios 
3.3  Influence  of  viscous  forces 

Scaling  hydraulic  models  is  dominated  by  Froude  number  scaling  (the  ratio  of 
gravitational  to  inertial  forces).  However,  there  are  other  important  forees  that  must  be 
considered. 

The  Reynolds  number  is  a measure  of  the  viscous  to  inertial  forces  in  a flow.  The  only 
way  to  maintain  the  effect  of  viscous  forces  is  to  ensure  that  both  the  Froude  number  and 
the  Reynolds  number  are  the  same  in  the  model  as  in  the  prototype.  This  necessarily 
involves  adjusting  fluid  viscosity  in  the  model  but  this  is  very  difficult  in  practice  and 
fresh  water  is  normally  the  only  fluid  suitable  for  modeling.  If  fluid  viscosity  is  the  same 
in  model  and  prototype  and  one  is  scaling  according  to  the  Froude  principles  it  becomes 
evident  that: 

Rer  = prVrLr/pr  = Lr'  ^ (22) 

This  means  that  Reynolds  numbers  will  be  much  smaller  in  the  model  than  in  the 
prototype. 

The  importance  of  maintaining  Reynolds  number  similarity  from  its  relationship  with 
another  important  dimensionless  variable,  the  drag  coeffleient  (Cd)  (equation  2)  which 
defines  the  relationship  between  drag  force  and  the  flow  dynamic  force. 

The  drag  coefficient  of  various  objects  is  shown  related  to  Reynolds  number  in  Figure  2. 
At  low  Reynolds  numbers  0.1<Re<1000  where  flow  is  laminar,  drag  on  bluff  bodies' 
(cylinders,  plates,  etc.)  increases  with  Re.  This  relationship  is  controlled  by  the  bodies’ 
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skin-friction,  that  is,  the  roughness  of  the  object's  surface  rather  than  its  overall 
dimensions.  At  higher  Reynolds  numbers  where  flow  is  turbulent,  skin  friction  ceases  to 
be  important,  and  it  is  the  overall  form  of  the  body  which  determines  the  drag  coefficient. 
Form  drag  remains  almost  constant  with  respect  to  Re  as  long  as  the  flow  is  turbulent.  It 
is  this  phenomenon  that  allows  the  scale  modeler  to  ‘cheat'  viscous  forces,  since  if  only 
form  drag  is  pertinent  in  the  experimental  situation,  the  drag  coefficient  will  remain 
roughly  similar  between  model  and  prototype  so  long  as  Reynolds  numbers  in  both  imply 
fully  turbulent  flow.  The  important  rule  in  modeling  is  that  Re  should  be  kept  as  high  as 
possible  in  the  model  situation,  but  need  not  match  Re  in  the  field.  Several  authors  give 
guidelines  on  the  minimum  Reynolds  number  required  for  accurate  modeling.  Allen 
(1947)  recommends  that  the  Re  in  the  model  should  not  be  less  that  1400. 


Figure  2:  The  relationship  between  Cd  and  Re  for  various  objects  (after  White, 
1994). 


The  Weber  Number  (see  equation  16)  defines  the  ratio  of  inertial  to  surface  tension 
forces  acting  on  a fluid.  Surface  tension  is  the  tensile  force  (N/m)  acting  perpendicular  to 
the  surfaee  of  a fluid  caused  by  intermolecular  forces.  Surface  tension  has  a negligible 
effect  upon  fluid  motion  in  most  rivers  because  it  is  a very  weak  force  compared  to  the 
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force  of  gravity.  However,  where  flow  depths  are  small,  as  in  the  case  of  hydraulic 
models,  surface  tension  may  distort  the  scaling  of  fluid  motion  characteristics  between 
prototype  and  model. 

There  is  no  consensus  to  suggest  a single  critical  Weber  number  for  models  to  avoid 
surface  tension  effects.  Studies  (see  Peakall  et  ah,  1996)  suggest  that  values  should 
exeeed  1 00.  White  ( 1 994)  suggests  that  the  Weber  number  is  important  only  if  it  is  of  the 
order  unity  or  less.  While  Novak  and  Cabelka  (1981)  suggest  three  minimum  operating 
conditions  to  minimise  surfaee  tension  effects.  These  are:  (I)  surface  waves  on  the 
model  should  be  greater  than  0.017m;  (2)  surface  flow  velocity  should  exceed  0.23  m/s  to 
allow  gravity  waves  to  develop,  and;  (3)  flow  depth  in  the  model  should  exceed  0.015  m. 
Noval  and  Cabelka  (1981)  also  state  that  the  minimum  Weber  number  should  be  11. 
Estimates  from  model  studies  range  from  2.5  to  17  and  from  29  to  160  (see  Peakall  and 
Warburton,  1996). 

3.4  Scaling  channel  boundary  resistance  in  fixed  bed  models 

In  the  previous  section  we  have  learned  the  importance  of  ‘Froude  scaling'  and  the 
scaling  laws  which  are  derived  as  a eonsequence.  We  now  need  to  apply  these  rules  to 
various  types  of  model  construction.  First  we  eonsider  models  where  one  is  simply 
interested  in  replieating  the  forees  and  flow  structure  of  a prototype.  Second  we  discuss 
the  far  more  complex  problem  of  scaling  deformable  boundary  models  and  similarity  of 
sediment  transport,  as  this  leads  us  towards  a solution  for  the  experimental  problem 
which  is  considered  in  this  researeh. 

It  is,  of  eourse,  evident  that  one  must  maintain  similarity  of  boundary  resistance  between 
model  and  prototype  ehannels.  Henderson  (1966)  uses  the  Manning  formula  to  express 
the  ratio  of  boundary  flow  resistance  between  model  and  prototype: 

vr= (23) 

nr 

where  Vr  = veloeity  ratio,  Rr  = hydraulic  radius  ratio,  Sr  = energy  gradient  ratio,  and  nr  = 
Manning's  n ratio. 
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If  we  intend  to  build  a model  which  has  geometric  similarity  it  follows  that  Sr=l  and 
Rr=Lr.  So,  rearranging  equation  23  with  respect  to  iir  and  bearing  in  mind  the  relation 
developed  in  equation  20  we  find  that; 


nr  - 


Lr 


1/2 


-=Lr 


1/6 


(24) 


1/6 


This  condition  must  be  met  if  dynamic  similarity  is  to  be  achieved.  The  relation  nr  = L 
is  very  similar  in  nature  to  Strickler's  definition  of  absolute  n which  he  derived  from  field 
evidence  in  gravel-bed  rivers  (Richards,  1982),  where: 


n = O.OI51D 


50 


1/6 


(25) 


This  similarity  lead  Einstein  and  Chien  (1954)  to  the  conclusion  that  n as  defined  by 
equation  25  for  a model  situation  is  directly  and  geometrically  proportional  to  n as 
empirically  measured  in  the  prototype.  In  other  words,  nr  is  dependent  upon  boundary 
geometric  similarity,  regardless  of  whether  friction  is  predominantly  caused  by  skin 
roughness  (as  will  be  the  case  in  the  model)  or  form  roughness  (as  will  be  the  case  in  the 
prototype).  We  can  therefore  combine  equations  24  and  25  to  give: 

nr^Dsor""  (26) 


Achieving  similarity  of  n is  not  quite  so  simple,  however,  because  n is  also  a function  of 
Reynolds  number  which  will  be  proportionally  much  smaller  in  the  model  as  compared  to 
the  prototype.  One  could  solve  this  problem  by  proportionally  reducing  surface 
roughness  in  the  model,  but  this  is  often  physically  impossible  if  the  prototype  surface  in 
question  is  already  smooth.  The  conclusion  is  that  turbulence  must  be  made  as  high  as 
possible  in  the  model  situation  in  order  to  preserve  boundary  roughness  scaling.  And  it  is 
usually  impossible  to  change  model  fluid  viscosity,  we  are  left  with  only  one  other 
variable  in  the  Reynolds  number  equation  which  we  can  adjust  to  keep  Re  as  high  as 
possible;  L or  flow  depth. 

3.5  Scale  distortion 
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Why  distort  model  geometric  scaling?  For  three  important  reasons:  ( 1 ) to  solve  the 
problem  of  low  model  Reynolds  numbers,  as  greater  flow  depths  in  a model  will  ensure 
that  the  flow  is  fully  turbulent;  (2)  because  pure  geometric  scaling  may  result  in  model 
flow  depths  that  are  too  small  to  enable  accurate  instrumentation,  and;  (3)  because  if 
depth  is  too  small  surface  tension  effects  may  distort  the  scaling  of  flow  processes 
between  model  and  prototype. 


If  a model  is  distorted,  Lr  must  be  split  into  two  variables;  Yr,  a vertical  scale,  and  Xr,  a 
horizontal  scale.  I'he  Froude  scaling  relationship  now  takes  the  form  (Henderson,  1966): 
Vr  = Yr"^  (27) 


because  it  is  vertical  rather  than  horizontal  distance  which  determines  the  effect  of 
gravity  on  velocity.  However,  velocity  acts  in  essentially  a horizontal  direction,  so  the 
time  scale  becomes: 


Tr  = Xr/Vr  = Xr/Yr"‘ 

(28) 

discharge  scale  becomes: 

Qr  = VrXrYr=XrYr''' 

(29) 

slope  scale  becomes: 

Sr  = Yr/Xr 

(30) 

and  force  scale  becomes: 

Fr=PrXr'YrYr/XrTr-'  = PrXr'Yr 

(31) 

Also,  substituting  equation  30  into  equation  23,  we  find  that: 

Rr2/3Y^l/2 


Which,  when  combined  with  equation  27  yields: 


nr  = 


Xr 


1/2 


(33) 


A consequence  of  equation  33  is  that  nr  will  usually  be  less  than  unity,  so  that  n in  the 
model  is  by  necessity  of  the  scale  distortion  greater  than  that  in  the  prototype,  which 
therefore  reduces  the  effect  of  low  Reynolds  numbers  in  the  model.  Correct  model  n 
values  are  usually  determined  on  a trial  and  error  basis  by  adjusting  roughness  until  the 
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stage-discharge  relationship  in  the  model  is  similar  to  that  of  the  prototype.  It  is  also 
worth  noting  that  in  some  instances  Sr  is  made  greater  than  Yr/Xr  in  moveable-bed 
models  in  order  to  assist  sediment  movement;  this  practice  is  known  as  'tilting'  the  model 
(Henderson,  1966). 

3.6  Moveable  bed  models 

Moveable  bed  models  are  more  complex  to  scale  than  fixed  bed  models  because 
parameters  such  as  sediment  density  and  grain  size  must  be  considered.  The  scale 
relationships  discussed  in  the  previous  section  must  also  be  maintained.  Moveable  bed 
models  were  built  on  a trial  and  error  basis  at  one  time  (Henderson,  1966).  Explicit 
design  methods  currently  available  (Yalin,  1971;  Einstein  and  Chien,  1954;  Blench, 
1955)  use  the  Shields  and  Particle  Reynolds  parameters  as  a basis  to  model  sediment 
transport.  The  rationale  for  this  is  that  the  Shields  and  Particle  Reynolds  numbers  are 
dimensionless  parameters  that  can  be  utilized  for  scaling  in  the  same  manner  as  the 
Froude  and  Reynolds  numbers. 

The  method  proposed  by  Einstien  and  Chien  (1954)  is  used  herein.  The  method  assumes 
that  the  shape  of  the  bed,  and  sediment  transport  properties  can  be  preserved  between 
model  and  prototype  by  exactly  matching  parameters  of  the  Shields  entrainment  function 
over  the  range  of  design  discharges.  Figure  3 shows  the  Shields  entrainment  function. 
The  Shields  function  is  defined  by  the  Shields  parameter  (0*),  where: 


6 


(34) 


where,  Ys  = bulk  unit  weight  of  sediment  = ps  (sediment  density,  kg/m^)  x g,  y = bulk  unit 
weight  of  water  = p x g (approx.  9810  N/m^),  D = representative  grain  size  (usually  taken 
to  be  Dso),  and  to  = boundary  shear  stress  , defined  as; 


To  = yRS 

and  the  Particle  Reynolds  number  (Re*),  where: 
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where  U*  = shear  velocity  = V xo  /p.  The  aim  is  to  make  0+p  / 0*m=l  and  Re*p/Re*m=l  in 
order  to  achieve  similarity  of  sediment  transport. 


Particle  Reynolds  Number  (R^,.) 

To  find  Tc  (critical  shear  stress)  calculate  D/v'/0,UYs/y-I  )gD  for  the  case  of  interest,  draw  a line 
corresponding  to  the  result  parallel  to  the  inclined  lines  and  read  off  the  value  of  0.  from  which 
Tc  can  be  calculated. 


Figure  3:  The  Shields  entrainment  function  (modified  from  Vanoni,  1966) 

If  bed  configuration  is  the  same  in  the  model  as  in  the  prototype,  boundary  roughness  (k) 
will  also  be  in  the  same  ratio  (k  = 1/30  x D65,  where,  D65  = Nikuradse  sand  roughness; 
Richards,  1982).  It  therefore  follows  that  n will  be  in  the  same  ratio,  so  we  can  substitute 
equation  26  into  33  to  give: 

R 

D30r"‘'=^  (37) 

that  solves  Manning  and  Froude  scaling  for  a mobile  boundary  scale  distorted  model. 


A second  consequence  of  boundary  roughness  being  in  the  same  ratio  as  model  and 
prototype  is  that  the  ratio  of  grain  resistance  (R’)  to  form  resistance  (R”)  will  equal 
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unity.  Thus,  boundary  shear  stress  will  also  be  proportional  between  model  and 
prototype,  so: 


TOr  = YrRrSr  = YrRrYr/Xr  (38) 

The  development  of  equations  37  and  38  mean  that  equations  34  and  36  can  now  be 
derived  in  terms  of  scale  ratios.  Thus  to  scale  sediment  transport  correctly  we  require: 
R.Y, 


(39) 


and: 


Re*r  - 


RrYrDr^ 

XrV'r- 


(40) 


Assuming  that  Vr  = I between  model  and  prototype  (fresh  water  in  both),  means  that  there 
are  four  independent  variables  (Yr,  Xr,  Dr,  and  ar)  in  equations  37,  39  and  40,  which  have 
to  be  solved  simultaneously  in  order  to  achieve  Froude  scaling  for  a mobile  boundary 
model.  The  designer  is  therefore  free  to  choose  one  of  the  four  independent  variable 
scale  ratios.  It  is  also  worth  noting  that  the  model  must  have  scale  distortion  (i.e.  Xr  ^ 
Yr)  in  order  to  solve  this  set  of  equations.  This  system  offers  an  excellent  analytical 
solution  to  mobile  boundary  scaling,  but  it  is  somewhat  inflexible  and  so  requires  that  the 
designer  have  unlimited  flume  construction  options  and  specialized  physical  materials  at 
his  disposal.  Scaling  was  pursued  in  the  flume  model  using  the  Manning-Froude  equation 
(37)  and  the  Shields-Particle  Reynolds  equations  (39  and  40). 
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4 Abiaca  Creek  Scaling  Method 

The  first  step  in  creating  a physical  model  was  to  select  a representative  reach  from  the 
study  sites  as  a prototype  for  the  model.  The  prototype  reach  selected  was  in  Abiaca 
Creek  in  Northern  Mississippi. 

The  appropriate  channel  dimensions  were  then  determined  by  assuming  that  the  most 
important  discharge  in  the  prototype  river,  in  terms  of  affecting  channel  geomorphology 
(the  dominant  discharge),  was  the  two-year  flood  (Qi).  Q2  is  often  considered  to  be  the 
dominant  geomorphic  discharge  (Summertleld,  1991),  and  the  rationale  for  this 
assumption  is  discussed  in  many  tluvial  geomorphology  texts  (Richards,  1982,  Knighton, 
1984). 


The  necessary  prototype  dimensions  were  obtained  from  survey  studies  conducted  by 
Watson  et.  al.  (1993)  and  are  listed  in  Table  3. 


Variable 

Symbol 

Value 

Two  Year  Discharge 

Q2 

48.1  nrVs 

Average  Width 

w 

17.9  m 

Average  Depth 

d 

1.9  m 

Hydraulic  Radius 

R 

1.57  m 

Width/Depth  Ratio 

w/d 

10 

Average  Bed-Slope 

S 

0.001 1 

Estimated  Manning's  n 

n 

0.032 

Sediment  D50 

D50 

0.25  mm 

Table  3:  Hydraulic  variables  for  the  prototype  reach  in  Abiaca  Creek  at  Q2 
discharge 

The  flume  used  to  conduct  the  experimental  study  had  an  overall  channel  width  of  0.63 
m,  a depth  of  0.61  m,  and  a total  length  of  10.0  m.  Slope  on  this  apparatus  was 
adjustable  by  means  of  a power-driven  hydraulic  jack.  Flow  entrance  conditions  were 
calmed  by  passing  the  water  vertically  through  a gravel  filled  feeder  tank  into  the  flume. 


As  discussed  above,  the  model  had  to  have  mobile  bed  and  banks  and  be  fully  Froude 
scaled.  However  it  was  also  neeessary  for  the  model  channel  boundary  to  be  at  the 
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threshold  of  motion  at  the  Q2  discharge  so  that  geomorphic  change  only  occurred  with 
the  introduction  of  simulated  LWD  elements.  In  this  way,  the  impact  of  LWD  on 
channel  geomorphology  could  therefore  be  isolated  from  other  channel  adjustment 
process.  The  prototype  channel  boundary  is  well  above  the  threshold  of  motion  at  the  Q2 
discharge  however,  so  it  was  not  possible  or  desirable  to  achieve  scaling  similarity  for 
sediment  transport  between  the  model  and  prototype.  Scaling  proceeded  using  the 
methodology  developed  by  Einstein  and  Chen  1954)  as  discussed  above. 

Froude  and  Manning  scaling  is  achieved  be  selecting  only  one  of  the  possible  hydraulic 
variables.  Given  that  the  flume  had  a width  constraint  of  0.61  m it  was  decided  to  fl.x  the 
initial  flume  width  at  0.3  m to  allow  for  lateral  adjustment  during  test  runs.  All  other 
model  scale  relationships  had  now  to  be  developed  to  solve  equation  39,  the  Froude- 
Manning  scaling  law.  We  therefore  required  Xri 
Xr  = Xp/X,r,=  17.9/0.3  = 59.667 
and  Yr.  If  the  model  had  no  scale  distortion  Yr  = Xr,  so: 

Ym  = Yp  / Yr  = 1 .9  / 59.667  = 0.03  m 

The  flume  depth  would  therefore  have  been  0.03  m.  However,  it  was  considered  that  this 
depth  was  too  small  as  fully  turbulent  flow  might  not  have  developed  and  surface  tension 
effects  could  well  have  been  significant.  The  solution  was  therefore  to  give  up  true 
geometric  scaling  by  distorting  the  vertical  scale.  This  practice  is  quite  legitimate  and 
was  an  absolute  necessity  for  equations,  37,  39  and  40  to  be  solved  simultaneously 
(Henderson,  1966).  We  doubled  the  geometrically  scaled  flume  depth  to  0.07  m.  Thus: 

Yr  = Yp/ Ym=  1.9/0.07  = 27.142 

Fastly  we  require  Rr,  which  must  be  derived  from  Yr  and  Xr.  Hydraulic  radius  (R)  is 
defined  as: 

R = A/P  (41) 

where  A = w x d,  and  P = wetted  perimeter  = 2d  + w.  So  the  hydraulic  radius  ratio  is 
determined  as: 

Rr  = YrXXr/[(2xYr)  + Xr]  (42) 

Rr  = 27. 1 42  X 59.667  / [(2  x 27. 1 42)  + 59.667]  =14.21 
Model  hydraulic  radius  is  therefore: 


L /.'.jM',;  '>■  .M'j  ' ■'  ''i‘ 

»(  '.n  I.' 


f 


V . « 


"1 


# 


■) 


I 


I 


i ,i 


I 


T 


I 


Rm  = Rp/Rr=  1.57/  14.21  =0.110m 
Rr  and  Xr  can  now  be  substituted  into  equation  22  to  yield  Dr: 

Dr'^^=  14.21^^^ / 59.667'^'  = 0.192 
The  required  boundary  grain  size  in  the  model  was  therefore: 

Dm  = Dp  / Dr  = 0.00025  / 0. 1 92  =1.3  mm 

Material  constraints  at  the  laboratory  however  meant  that  the  closest  grain  size  available 
was  0.8-mm  quartz  sand,  which  had  a density  of  2620  kg/m  . The  model  boundary 
roughness  was  therefore  slightly  lower  than  that  required  for  pure  Froude-Manning 
scaling,  the  consequence  being  that  the  model  channel  had  slightly  less  boundary 
resistance  than  the  prototype. 

Now  that  Xr  and  Yr  were  calculated  it  was  possible  to  determine  the  required  tlume  slope, 
Sr: 

Sr  = Yr  / Xr  = 27. 1 42/  59.667  = 0.45 

Thus: 

Sm  = Sp  / Sr  = 0.00 1 1 / 0.45  = 0.0024  m/m 

At  this  point  the  model  dimensions  were  scaled  as  accurately  as  physieal  eonstraints 
would  allow  aceording  to  Froude-Manning  scaling  laws  for  a depth-distorted  model. 

Given  these  initial  dimensions,  tests  were  conducted  to  determine  what  flume  geometry 
was  practieally  possible,  given  that  the  channel  boundary  was  to  be  molded  from  0.8-mm 
sand.  Through  trial  and  error  tests  it  was  found  that  the  stable  angle  of  repose  for  the 
saturated  sand  was  34°  degrees.  The  channel  banks  could  therefore  not  be  made  near 
vertieal,  as  was  the  case  in  the  prototype.  The  elosest  approximation  to  the  prototype  for 
the  flume  geometry  as  required  by  the  scaling  parameters  that  had  stable  banks,  is  shown 
in  Figure  4. 
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Figure  4:  Cross  section  of  flume  channel  (dimensions  in  meters). 

This  channel  has  the  following  hydraulie  dimensions:  (1)  cross-sectional  area  = 0.014  m^, 
(2)  wetted  perimeter  = 0.35  m,  (3)  hydraulic  radius  = 0.04  m. 

It  should  be  noted  that  the  measured  hydraulic  radius  for  the  tlume  (0.04  m)  is  smaller 
than  the  hydraulic  radius  determined  through  the  scaling  technique  (0.11  m).  This 
difference  is  due  to  the  fact  that  the  channel  had  to  have  sloping,  rather  than  vertical 
banks.  This  difference,  however,  does  not  significantly  compromise  the  initial  scaling 
procedure.  Indeed,  the  discrepancy  between  the  actual  and  scale-derived  flume  R even 
helped  to  maintain  parity  of  effective  boundary  roughness  between  flume  and  prototype 
by  compensating  for  the  mismatch  between  the  boundary  grain  size  required  through 
scaling  and  the  smaller  grain  size  which  had  to  be  used. 

The  model  was  not  complete  at  this  stage,  however,  because  we  required  that  the  channel 
boundary  be  at  the  threshold  of  motion.  Mobile  boundary  scaling  techniques  now  had  to 
be  employed.  Adopting  the  rational  scaling  method  proposed  by  Einstein  and  Chien 
(1954),  solutions  had  to  be  found  to  equations  41  and  42  for  the  model-prototype 
dimensions.  First,  scale  relationships  had  to  be  determined  for  the  fluid  kinematic 
viscosity  ratio  (vr)  and  sediment  density  ratio  (yr)  in  order  to  solve  these  equations. 

Given  that  fresh  water  had  to  be  used  in  the  model  Vr  was  assumed  to  equal  one: 

Vr  = Vp/Vm  = 1 X 10'^/  1 X 10'^=  1 
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sediment  density  of  bed  material  in  the  prototype  is  roughly  2620  kg/m^  (Watson  et  al„ 
1993),  which  is  the  same  as  that  of  the  sand  used  in  the  model,  so  ar  was  also  assumed  to 
equal  one: 

Yr  = Yp/Ym  = 2620  x 9.8l  /2620  x9.81  = l 

The  necessary  scale  ratios  were  all  derived  so  it  was  now  possible  to  solve  the  Shields 
and  Particle  Reynolds  scaling  functions.  All  scale  ratios  and  physical  dimensions  for  the 
flume  model  are  shown  in  Table  4. 


Scaling  factor 

Symbol 

Scaling 

value 

Comparable 

dimension 

Flume  value 

Horizontal 

Xr 

59.667 

Flow  width 

0.300  m 

Vertical 

Yr 

27.142 

Flow  depth 

0.070  m 

Hydraulic  radius 

Rr 

14.212 

Hydraulic  radius 

0.04  m 

Slope 

Sr 

0.45 

Water  surface  slope 

0.0024 

Kinematic 

Vr 

1 

Kinematic 

1x10'^  m“/s 

viscosity 
Bulk  sed.  density 

Yr 

1 

viscosity 
Bulk  sed.  weight 

25702  N/m^ 

Sediment  grain  size 

Dr 

0.192 

Sediment  D50 

1 .3  mm 

Table  4:  Scaling  factor  values  and  main  flume  dimensions  determined  by  Froude- 
Manning  sealing 


The  threshold  boundary  condition  desired  for  the  flume  model,  the  Q2  discharge,  was  not 
similar  to  the  boundary  condition  in  the  prototype  reach  at  Q2  because  as  significant 
sediment  transport  occurs  in  Abiaca  Creek  even  at  flows  less  than  Q2.  Equations  39  and 
40  could  not  therefore  have  products  of  unity  to  fulfil  the  requirements  of  the  model.  It  is 
instructive  however  to  calculate  the  Shields  parameter  and  Particle  Reynolds  parameter 
ratios  for  the  scale  ratios  presented  in  Table  4 for  demonstration. 


Remember  that  0*r  = RrYr/yrXrDr.  Thus: 

e*r  = (14.212  X 27.142)  / (lx  59.667  x 0.192)  = 33.667 
And  that  Re*r  ^ RrYrDr^/XrUr^.  Thus: 

Re*r  = (14.212  X 27.142  x 0.192^ )/  (59.667  x 1“  )=  0.238 
We  see  to  restore  0*r  to  unity  without  distorting  the  Rr,  Xr,  Yr  and  Dr  determined  by 
Froude  scaling,  yrimist  be  made  greater  than  one  by  reducing  the  model  sediment  density. 
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To  restore  Re*r  to  unity,  Vr  must  be  made  less  than  one  by  increasing  the  model  fluid 
viscosity.  These  adjustments  were  not  possible  given  the  laboratory  equipment  available. 

The  hydraulic  conditions  at  the  threshold  of  motion  must  be  determined.  First,  we  had  to 
determine  the  Shields  and  Particle  Reynolds  numbers  for  the  hydraulic  variables  defined 
in  Table  4,  and  find  where  the  two  values  intersected  on  the  Shields  entrainment  function 
diagram  (see  Figure  3).  These  parameters  were  calculated  as  follows  (note  that  we  have 
used  the  measured  channel  hydraulic  radius  of  0.04  m,  rather  than  the  theoretical  value 
determined  by  scaling  (0.1 1 m)  as  the  former  value  represents  the  true  model  hydraulic 
condition).  The  Shields  parameter  = 0*=To/(ys-y)D5o,  where,  to=yRS.  Thus: 

= 98 1 0 X 0.04  X 0.0024  / [(25702  - 98 1 0)  x 0.0008]  = 0.073 
and  the  Particle  Reynolds  parameter  = R^;*  = U^Dso  / v,  where,  U+  = Vxo  / Pw  Thus: 

Rc*m  = V((98 1 0 X 0.04  X 0.0024)  / 1 000)  x 0.0008  / 1 x 1 0'^  = 24.55 
Plotting  these  values  on  the  Shields  entrainment  function  (Figure  3),  the  intersect  is 
above  the  critical  curve  (square  marker)  so  the  model  did  not  have  a threshold  boundary 
condition  for  the  pure  Froude-Manning  scaled  variables.  We  therefore  employed  the 
technique  of  Tilting'  the  model  by  adjusting  the  prescribed  slope  so  that  the  Shields 
parameter  was  reduced  to  the  critical  value  (0+c)-  The  model  therefore  had  a slope 
condition  where  Sr  Yr  / Xr. 

The  slope  required  for  a threshold  boundary  condition  was  determined  by  substituting  an 
appropriate  Shields  value  into  equation  36  and  solving  for  the  slope.  The  threshold  curve 
in  Figure  3 can  be  used  to  visually  determine  a critical  Shields  value  for  a known  particle 
Reynolds  number.  However,  methods  developed  by  Bogardi,  (see  Henderson,  1966),  and 
Vanoni  (1966)  more  accurately  determine  0*c  for  a given  boundary  roughness  based  upon 
empirically  derived  equations.  In  Vanoni's  method  Q*c  is  determined  by  calculating  a 
value  for  the  following  equation: 


0.1  ^-1  g03, 
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A line  is  then  drawn  eorresponding  to  that  value  parallel  to  the  inelined  lines  on  Figure  3, 
down  to  the  threshold  curve,  and  then  horizontally,  parallel  to  the  x axis,  to  intersect  with 
0*c  on  the  y axis.  We  used  this  procedure  to  determine  0+c  for  the  model.  Substituting 
model  values  into  equation  43  gives: 

0.0008/1  X 1 0'S  a/O.  I (25702  / 98 1 0)  - 1 ) 9.8 1 x 0.0008  = 2.85 
Drawing  a line  corresponding  to  this  value  parallel  to  the  inclined  lines  on  Figure  3 and 
reading  off  the  value  horizontally  yielded  a critical  Shields  parameter  0*c  of  0.060.  This 
value  was  then  substituted  into  the  Shields  equation  to  determine  the  critical  slope  for  the 
model: 

0.060  = 98 1 0 X 0.04  x S^  / (25702  - 98 1 0)  x 0.0008 
which  yields: 

Sn,- 0.0019 

It  was  therefore  necessary  to  initially  relax  the  flume  slope  from  the  Fronde  scaled  value 
of  0.0024  down  to  0.0019  in  order  to  guarantee  a critical  threshold  boundary  condition. 

The  flume  channel  was  then  constructed  using  the  dimensions  shown  in  Figure  4,  with  a 
bed  slope  of  0.0019.  A discharge  was  then  run  through  the  channel  and  increased 
incrementally  until  flow  depth  reached  0.07  m.  The  flume  was  then  left  to  stabilize  for 
several  hours  to  make  sure  that  sediment  transport  did  not  occur.  No  transport  was 
observed  so  slope  was  increased  slowly  to  determine  where  the  critical  condition  lay. 
Through  trial  and  error  it  was  discovered  that  the  water  surface  slope  required  for  the 
threshold  of  particle  motion  was  actually  0.0022,  a value  close  to  the  slope  of  0.0024  that 
was  initially  required  for  scaling  similarity.  Consequently,  the  channel  slope  was  set  to 
0.0022  giving  a pre-test  channel  that  closely  achieved  both  Froude-Manning  scaling  and 
a threshold  boundary  condition. 

At  the  bankfull  stage,  with  channel  slope  set  to  0.0022,  model  discharge  (Qm)  measured 
at  the  flume  outlet  was  0.0033  m^/s. 

This  measured  value  can  be  compared  with  the  model  discharge  predicted  by  the  scaling 
laws  (Oin(s))  to  provide  a check  on  the  calculations.  Remember  that  Qr  = XrYr'  ^ (equation 
29).  Now  this  equation  assumes  continuity  of  flow  in  the  manner: 
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Q = V X A 


(44) 


where  A = wxd  = XxY. 

However,  the  channel  cross  section  in  the  model  is  not  rectangular,  so  A X x Y.  Cross- 
sectional  area  in  the  tlume  is  0.014  m,  so  it  is  a simple  matter  to  calculate  a critical  depth 
(Yc)  to  satisfy  equation  44. 

Yc  = A / X = 0.014  / 0.3  = 0.046  m 
Now  calculating  Yr  using  Yc  we  find: 

Yr=  1.9/0.046  = 41.30. 

Qr  can  now  be  calculated  correctly: 

Qr=  59.667  x41.30'^=  15836.5 
Model  discharge  as  determined  by  scaling  (Qm(s))  is  therefore: 

Qm(s)  = Qp  / Qr  = 48. 1 / 1 5836.5  = 0.0030  m^/s 
Discharge  derived  by  scaling  (0.0030  mVs)  was  therefore  in  good  agreement  with  the 
discharge  measured  at  the  flume  outlet  (0.0033  m'^/s),  and  we  were  therefore  satisfied  that 
all  variables  had  been  calculated  accurately. 

Finally,  depth  averaged  velocity  was  derived  for  the  flume  (Vm)  using  equation  44.  This 
yields: 

Vm  = 0 / A = 0.0033  / 0.0 1 4 = 0.24  m/s. 

A comparison  of  Froude  values  in  the  prototype  channel  and  model  channel  should  yield 
the  following  relationship  if  scaling  is  correct:  Frp  / Frm  = 1 (equation  20).  Using  the 
derived  flow  velocity  of  0.24  m/s  the  Froude  number  for  the  model  was: 

Frm  = V / Vgd  = 0.24  / V9.8 1 x 0.07  = 0.29 
While  the  Froude  number  for  the  2-year  discharge  event  in  Abiaca  Creek  is: 

Frp=  1.4/ V9.81  X 1.9  = 0.32 

Thus: 

Frp  / Frm  = 0.32  / 0.29  = 1.10 

The  model  flow  conditions  were  therefore  very  close  to  the  ideal  condition  for  unity  in 
Froude  scaling. 
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Reynolds  number  for  the  tlume  was: 

Rem  = vd  / V = 0.24  0.07  / 1 x 1 0'^  = 1 6800 

while  Reynolds  number  for  the  2-year  discharge  event  in  Abiaca  Creek  is: 

RCp  = 1 .4  X 1 .9  / 1 X 1 0'^  = 2660000 

There  is  therefore  a great  difference  in  the  model  and  prototype  Reynolds  numbers. 
However,  because  the  model  Reynolds  number  is  high  enough  to  ensure  fully  turbulent 
flow  (Allen  (1947)  suggests  model  values  should  e.xceed  1400),  the  difference  in  absolute 
values  is  not  critical.  This  is  because  drag  coefficients  on  immersed  bodies  remain 
almost  identical  for  the  entire  range  of  flow  Reynolds  numbers  that  signify  a fully 
turbulent  flow  (refer  to  Figure  2).  The  flume  model  was  considered  adequately  scaled 
according  to  the  two  most  important  hydraulic  dimensionless  parameters. 

One  further  dimensionless  coefficient,  the  Weber  number,  had  also  to  be  determined  for 
the  flume  however  in  order  to  ensure  that  surface  tension  effects  were  not  significant. 
Recall  that  the  Weber  number  = We  = v“  p^d/  cr.  Thus: 

Wcm  = 0.24“  X 1 000  X 0.07  / 0.073  = 55.23 

This  value  lie  in  the  middle  of  the  range  quoted  for  other  hydraulic  models,  while  the 
criteria  suggested  by  Novak  and  Cabelka  ( 1981 ) are  all  exceeded  by  the  model  (although 
length  of  surface  waves  were  unknown).  So,  despite  some  confusion  in  the  literature 
concerning  a precise  limiting  Weber  number,  we  were  satisfied  that  the  model  at  least 
met  the  minimum  criteria  to  ensure  that  fluid  motion  would  was  not  significantly  affected 
by  surface  tension  effects  during  the  test  runs.  Table  5 shows  a comparison  of  prototype 
and  model  variable  dimensions.  Figure  5 shows  a flow  diagram  of  the  scaling  procedure 
used  to  construct  the  model. 
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Variable 

Symbol 

Units 

Model 

Prototype 

Channel  width 

w 

m 

0.30 

17.90 

Reach  length 

L 

m 

1.00 

65.00 

Flow  depth 

d 

m 

0.07 

1.90 

Bed  slope 

S 

0.0022 

0.001 1 

Hydraulic  radius 

R 

m 

0.04 

1.57 

Cross-sectional  area 

A 

m^ 

0.014 

34.00 

Wetted  perimeter 

P 

m 

0.35 

21.70 

Bankful  1 discharge 

Qb/Q2 

m^/s 

0.0033 

48.10 

Derived  mean  flow  velocity  at 

V 

m/s 

0.24 

1.40 

Qb 

Representative  grain  size 

D50 

mm 

0.8 

0.25 

Sediment  density 

Ps 

kg/m^ 

2650.0 

2620.0 

Width-to-depth  ratio 

w/d 

- 

4.20 

10.00 

Particle  Reynolds  number  at  Qb 

R,* 

- 

2.32 

3.25 

Shields  Parameter  at  Qb 

0 + 

- 

0.060 

1.27 

Shear  velocity  at  Qb 

U* 

- 

0.029 

0.130 

Manning's  n 

n 

- 

0.023 

0.032 

Important  Dimensionless  Coefficients 

Froude  number 

Fr 

- 

0.29 

0.32 

Reynolds  number 

Re 

- 

16800 

2660000 

Weber  number 

We 

- 

55.23 

51013.70 

Table  5:  Comparison  of  model  and  prototype  channel  variables 
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Figure  5:  Flow  diagram  of  procedure  for  creating  a Froude  scaled  mobile  boundary 
flume  model  from  a prototype  river 
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Figure  5:  Continued 
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5 Flume  Design  and  Construction 

The  intention  of  the  flume  analysis  was  to  test  the  geomorphic  and  hydraulic  impact  of 
the  four  distinct  jam  types  (Underflow,  Dam,  Deflector,  and  Flow  Parallel)  which  were 
identified  in  the  field  and  conceptualized  in  the  Debris  Jam  Classification  Model  (refer  to 
Wallerstein  et  al.,  1997).  Two  important  simplifications  were  made  in  the  flume  model 
to  ensure  that  results  would  be  repeatable.  First,  rather  than  creating  jams  made  from  a 
number  of  complex  elements,  single  cylindrical  elements  were  used  representing  the 
‘key*  debris  element  in  each  jam  type.  Second,  rather  than  using  one  average 
representative  key  debris  (tree)  height  and  varying  the  flume  channel  width  to  simulate 
the  different  debris  length  to  channel  width  ratios  suggested  in  the  Debris  Jam 
Classification  Model,  the  reverse  was  done.  Flume  channel  width  was  fixed  and  debris 
element  lengths  were  changed  to  simulate  varying  debris  length  to  channel  width  ratios. 
This  compromise  was  made  to  maintain  scale  similarity  in  the  flume.  Table  6 details  the 
relative  and  absolute  proportions  of  the  LWD  elements  used  in  the  flume. 


Each  element  was  positioned  in  the  channel  to  simulate  the  natural  location  of  LWD  if  it 
were  to  fall  perpendicular  to  the  flow  in  the  prototype  channel.  Figure  6 shows  a 
schematie  diagram  of  the  element  positions  in  the  flume. 


Debris  Jam 

Classification 

Type 

Average  debris 
lengths  as  a 
proportion  of 
channel  width 

LWD  element 
lengths  in 
flume  (m) 

Average  LWD 
diameter  as  a 
proportion  of 
flow  depth 

LWD  element 
diameter  in 
flume  (m) 

Underflow 

1.00 

0.30 

0.27 

0.019 

Dam 

0.50 

0.15 

0.27 

0.019 

Deflector 

0.33 

0.10 

0.27 

0.019 

Parallel 

0.23 

0.07 

0.27 

0.019 

Table  6:  LWD  proportions 


The  0.3  element,  representing  underflow  type  LWD,  almost  spans  the  channel  at  the 
bankful  level.  The  0.15  m,  0.1  m and  0.07  m elements  were  positioned  to  simulate  the 
hypothetical  case  of  LWD  falling  from  the  left  bank  into  the  channel.  It  must  be  noted 
that  the  elements  were  not  touching  the  channel  boundary  at  the  start  of  each  run. 
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The  flume  apparatus  is  displayed  in  Figure  7.  The  channel  was  created  by  pre-wetting 
the  sand  and  then  cutting  the  cross-sectional  form  into  the  length  of  the  flume  using  a 
scraper  plate,  machined  to  the  dimensions  given  in  Figure  4,  and  attached  to  a carriage 
mounted  on  top  of  the  flume.  Excess  sand  was  carefully  removed  and  any  irregularities 
in  the  channel  carefully  smoothed.  Plate  1 shows  the  model  channel  and  scraper  plate 
prior  to  running. 

Slope  was  adjusted  to  the  correct  level  using  a calibrated  Jack  attached  to  the  upstream 
end  of  the  flume.  An  adjustable  recirculating  pump  was  used  to  drive  the  flow  through 
the  flume,  and  discharge  was  calculated  using  a v-notch  weir.  The  weir  was  set  into  the 
outlet  of  a collecting  tank  located  beneath  the  flume  outlet. 

It  was  necessary  to  install  a rectangular  weir  plate  at  the  flume  outlet  which  could  be 
raised  from  the  bed  upwards  in  order  to  prevent  the  flow,  accelerating  over  the  outlet, 
from  increasing  the  total  head  in  the  lower  end  of  the  tlume.  This  discontinuity  in  the 
energy  gradient  caused  limited  bed  load  sediment  transport  towards  the  end  of  the  flume 
until  the  problem  was  rectified  using  the  weir  plate.  Plate  2 shows  the  undisturbed  water 
surface  in  the  channel,  upstream  of  the  test  section. 

Debris  elements  were  made  from  cylindrical  aluminum  stock,  cut  to  appropriate  lengths. 
A mechanical  strain  gauge  was  devised  to  calculate  the  dynamic  pressure  force  on  the 
elements  during  test  runs.  This  apparatus  is  shown  in  Figure  8. 
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0.07m  Element 

Planform  Cross  Section 


0.1  m Element 

Phmform  Cross  Section 


0.15  m Element 

Planform  Cross  Section 


0.3  m Element 


Planform  Cross  Section 


0 1 10 
Units  are  in  centimeters 


Eigure  6:  Schematic  diagrams  of  LWD  positions  in  the  flume 
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Instrument  carriage 
mounted  on  larger  flume 


Scraper  plate  to  Balance  / 
cut  channel 


Pivot  point  fixed 
to  flume  top 


Entrance  plate 


Flow  inlet 
Adjustable  bed  slope 


0.8  mm  sand 


Exit  plate  ^ 

Pivot  point  on  outer  flume 


Figure  7:  Flume  set-up 


EWD  element 


Figure  8:  Strain  gauge  set-up 
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Each  LWD  element  had  a thread  drilled  at  the  center  of  its  long  axis  which  enabled  it  to 
be  screwed  into  a steel  spindle  which  was  in  turn  attached  to  a swinging  arm  that  was 
suspended  by  knife  edge  pivots  in  a eradle  attaehed  to  the  flume  top.  Element  height 
within  the  flow  could  be  adjusted  by  sliding  the  spindle  up  and  down  within  the  pivoting 
arm  between  two  locking  nuts. 

If  the  arm  and  spindle  were  free  moving,  force  on  the  element  would  push  it  forward  in 
the  flow,  the  arm  rotating  about  the  pivot  point.  However,  the  pressure  foree  was 
translated  to  a strain  line  made  from  fine  fishing  wire,  whieh  eonneeted  to  a loading  hook 
in  the  base  of  an  electronic  balance.  Thus,  foree  on  the  element  was  translated  to  a mass 
reading  on  the  balance.  An  equivalent  foree  reading  was  caleulated  using  the  equation; 

F|=G|/kgxg  (45) 

where  Fi  = equivalent  force  reading  on  the  balanee  (N),  Gj  = mass  reading  on  the  balance 
(grams),  kg  = conversion  to  kilograms  (divide  by  1000).  Actual  drag  force  on  the  element 
was  determined  by  calculating  the  moment  of  force  about  the  pivot  point.  Thus: 

F2  = (F, -Fo)  X, /(X, +X2)  (46) 

where  F2  = drag  force  on  FWD  element  (N),  Fo  = residual  force  reading  on  balanee  when 
there  is  no  flow  due  to  pre-loading  tension  in  the  strain  line  (N),  X|  = distance  from  pivot 
point  to  strain  line  (m),  and  X2  = distance  from  strain  line  to  center  of  element  (m). 

For  each  test  run  the  channel  was  re-eut,  smoothed  and  the  design  discharge  run  through 
for  two  hours  prior  to  testing  to  ensure  that  the  boundary  remained  stable.  The  element 
and  strain  gauge  apparatus  were  then  set  up  on  the  flume  and  the  pivot  arm  lowered  into 
its  cradle,  so  that  the  FWD  element  became  suspended  in  the  flow  at  the  appropriate 
depth  and  position  relative  to  the  ehannel  eross-section.  Plate  3 shows  the  strain  gauge 
apparatus  mounted  on  the  flume  during  a test  run. 
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Plate  3:  Experimental  set-up,  showing  pivot  arm  and  balance.  Run  with  0.1  m 
element 


Readings  on  the  balance  were  then  taken  as  soon  as  possible  after  element  immersion  and 
then  at  fixed  time  intervals  which  lengthened  gradually  as  force  readings  declined 
asymptotically  due  to  channel  re-adjustment  to  an  equilibrium  form.  Adjustment  of  the 
channel  geomorphology  was  monitored  during  each  run  using  a video  suspended  above 
the  flume  section  and  also  noted  by  sight. 

The  largest  element  caused  erosion  and  sediment  transport  that  finally  ceased  after  four 
hours  of  running  time.  Force  readings  therefore  also  decayed  to  a roughly  constant  value 
after  this  time  period.  The  remaining  elements  were  tested  in  order  of  size.  Time  until 
sediment  transport  ceased  was  found  to  decrease  with  element  size  but  total  monitoring 
time  was  maintained  at  four  hours  to  ensure  integrity  in  the  resulting  data-set.  At  the  end 
of  each  test  run  the  flume  was  drained  and  a detailed  survey  made  of  the  channel 
topography  in  a section  1.2  meters  long  downstream  from  the  element  suspension 
position,  using  a point  gauge  mounted  on  the  flume  instrument  carriage. 
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6 Description  of  the  Flume  Results 

6.1  Introduction 

There  were  two  principle  aims  in  the  flume  experiments:  (1)  to  determine  the  scale,  type 
and  spatial  extent  of  adjustment  in  channel  geomorphology  caused  by  Underflow,  Dam, 
Deflector  and  Flow-Parallel  debris  elements  in  isolation  from  other  channel  adjustment 
processes  such  as  degradation,  and  (2)  to  understand  the  hydraulic  processes  associated 
with  simple  LWD  forms.  Quantification  of  the  hydraulic  processes  provides  insight  into 
the  affect  of  LWD  on  flow  field,  velocity  and  shear  stress  within  a channel  reach,  which 
in  turn  drive  erosion  and  sediment  transport  processes  thus  linking  LWD  with  specific 
geomorphic  forms. 

The  forces  acting  on  LWD  elements  can  be  determined  by  re-arranging  the  dimensionless 
drag  coefficient  equation  (equation  2)  to  give: 

Fd  = 0.5  v^  Ad  Pw  Cd  (47) 

where  Fd  = Drag  force  (N). 

There  are  three  unknown  quantities  in  the  equation  for  a given  LWD  element  (Ad  is 
known).  Force  was  measured  directly  in  the  experiments  thus  making  it  necessary  to 
determine  a drag  coefficient  for  each  element  in  order  to  solve  the  equation  for  approach 
flow  velocity. 

6.2  Expansion  of  flume  results  to  prototype  dimensions 

The  key  to  these  flume  experiments  is  that  the  model  is  scaled  from  a prototype  river 
reach,  maintaining  dynamic,  kinematic,  and  geometric  similarity.  It  is  therefore  a 
legitimate  procedure  to  re-scale  the  model  results  using  the  various  scaling  factors 
described  earlier  to  produce  force,  velocity,  time  and  geomorphic  change  dimensions  for 
the  prototype  reach.  By  doing  so,  the  hydraulic  and  geomorphic  impact  of  LWD  on  sand 
bed  rivers  can  be  described. 
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The  variables  that  required  scaling  to  prototype  dimensions  from  the  model  were:  (1) 
time  (T,t,  ^ Tp),  (2)  velocity  (Vm  ^ Vp),  and  (3)  drag  Force  (Fm  ^ Fp).  Time  data  is 
scaled  according  to  equation  28.  Time  scale  ratio  (Tr)  is  therefore: 

Tr  = (59.667/27. 142)°^=  1.48 

Time  in  the  prototype  (Tp)  is  thus  derived  from  time  in  the  model  (Tm)  by: 

Tp  = Tn,x  1.48 

Velocity  data  is  scaled  according  to  equation  27.  The  velocity  ratio  is  therefore: 

Vr  = 27.142°-^  = 5.2 

Velocity  in  the  prototype  (Vp)  is  therefore  derived  from  velocity  in  the  model  (Vm)  by: 

Vp  = Vm  X 5.2 

Drag  force  is  scaled  according  to  equation  3 1 . Water  density  in  the  model  and  prototype 
were  assumed  to  be  similar  so  that  the  water  density  ratio  (pr)  is  equal  to  one.  Drag  force 
ratio  (Fr)  is  therefore: 

Fr  = 1 X 59.667^  X 27.142  = 96629.6 

Force  in  the  prototype  (Fp)  is  therefore  derived  from  force  data  in  the  model  (Fm)  by: 

Fp  = Fm  X 96629.6 

These  scaling  values  were  then  used  to  convert  model  results  to  the  prototype 
dimensions. 

The  flume  model  debris  element  dimensions  also  have  to  be  re-scaled  to  prototype 
dimensions  in  order  to  make  sense  of  the  model  data.  Using  the  horizontal  and  vertical 
scaling  relationships  developed  in  section  4 dimensions  for  the  prototype  debris  elements 
are  listed  in  Table  7. 

6.3  Limitations  of  the  flume  model 

There  are  three  important  differences  between  the  flume  model  and  the  prototype  that 
must  be  considered  when  interpreting  the  data.  The  prototype  has  a considerable 
sediment  transport  load  through  the  channel  at  the  Q2  discharge,  whereas  the  flume 
model  had  no  sediment  transport,  except  for  that  caused  by  the  LWD  elements.  In  the 
field  this  transported  sediment  will  obviously  interact  with  LWD  and  may  be  deposited  as 
bar  formations  in  areas  of  the  channel  where  the  debris  shields  the  flow,  or  causes 


41 


p 


I ) , I . N 


I 


''r'-'J.-g 


backwater  effects.  Such  formations  were  observed  in  the  field.  As  there  was  no 
sediment  transport  from  upstream  in  the  flume  it  is  therefore  likely  that  the  volume  and 
areal  extent  of  bars  created  by  each  element  would  be  less  than  those  observed  in  the 
prototype. 


Representative 
Element  Type 

Element 
Length  in 
Model  (m) 

Element 
Length  in 
Prototype 
(m) 

Element 
Diameter 
in  Model 
(m) 

Element 
Diameter  in 
Prototype 
(m) 

Underflow 

0.30 

18.0 

0.019 

1.1 

Dam 

0.15 

9.0 

0.019 

1.1 

Deflector 

0.10 

6.0 

0.019 

1.1 

Flow-parallel 

0.07 

4.2 

0.019 

l.l 

Table  7:  Prototype  debris  element  dimensions 


The  composition  of  the  banks  of  the  prototype  river  is  a mixture  of  materials  with 
different  structural  properties.  The  banks  of  the  base  flow  channel  are  often  found  to  be 
non-cohesive  sand,  if  the  reach  is  stable  or  aggrading,  or  highly  cohesive  clay  material  if 
the  channel  is  degrading.  Bank  angles  in  this  region  tend  to  be  shallow  (<40°).  These 
materials  are  fluvially  eroded  through  particle-by-particle  entrainment.  Above  the  base 
flow  area,  banks  are  often  found  to  be  composed  of  loess  or  a variety  of  different  alluvial 
deposits.  Bank  angles  in  these  materials  tend  to  be  close  to  vertical.  These  higher  areas 
of  the  channel  bank  come  into  contact  with  the  flow  only  during  the  less  frequent  events 
(remember  if  the  creek  has  degraded,  sueh  events  are  unlikely  to  be  ‘bankfulf),  and  tend 
to  retreat  mainly  by  mass  failure  as  the  lower  bank  material  is  fluvially  eroded  away. 
LWD  tends  to  cause  an  increase  in  fluvial  erosion  on  the  bed  and  of  the  lower  banks. 
This  in  turn  results  in  mass  failure  of  the  upper  bank,  so  that  large  block  of  material  fall 
into  the  ehannel,  onto  and  around,  the  debris  aceumulation.  The  flume,  however,  had 
banks  composed  solely  of  sand,  with  an  consistent  angle  of  34  degrees.  As  accurate 
representation  of  more  cohesive  materials  at  that  scale  is  virtually  impossible.  So,  while 
the  flume  banks  were  of  an  angle  and  consistency  roughly  similar  to  that  found  in  the  low 
flow  region  of  stable  or  aggrading  channels,  they  were  not  representative  of  degrading 
channel  environments,  or  bank  material  properties  generally  found  higher  in  the  bank 
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profile.  LWD  induced  bank  erosion  in  the  tlume  therefore  only  occurred  through  direct 
fiuvial  action,  as  it  was  not  possible  to  replicate  the  mass  failure  erosion  mechanism.  The 
resulting  channel  morphologies  created  in  the  flume  will  therefore  only  partially 
represent  the  forms  that  are  found  in  the  field,  although  we  consider  that  bed  features  and 
lower  bank  erosion  patterns  will  have  been  faithfully  reproduced. 

LWD  jams  in  the  field  are  complex.  They  have  generally  been  found  to  form  around  a 
'key'  debris  element  such  as  a large  tree.  They  also  evolve  over  time  becoming  either 
larger  and  more  impermeable  as  smaller  debris  transported  from  upstream  piles  up 
against  the  Jam.  The  'jams'  in  the  flume,  however,  were  very  basic  and  represented  'key' 
debris  element  sizes  and  positions  in  the  channel  for  the  generalized  Underflow,  Dam, 
Deflector,  and  Flow  Parallel  jam  classifications.  It  would  have  been  virtually  impossible 
to  measure  forces  on  a 'realistic'  jam  with  a ramification  of  smaller  and  more  complex 
element  forms.  Moreover,  the  consistency  of  element  form  (except  for  the  all  important 
length  dimension)  made  it  possible  to  attribute  the  difference  in  hydraulic  and 
geomorphic  results  solely  to  the  variation  in  area  of  flow  blocked  by  the  element. 

It  is  virtually  impossible  to  create  a scale  model  of  any  natural  channel  environment 
without  resorting  to  a considerable  level  of  simplification.  The  simplification  was 
necessary  as  the  overall  aim  of  these  experiments  was  to  understand  the  gross  properties 
of  river  flow  interaction  with  different  sizes  of  tree  like  objects  and  the  resulting  potential 
effects  had  upon  the  channel  boundary. 
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7 Discussion  of  Results 


7.1  Drag  coefficient  calculations  for  simulated  debris  elements 

Drag  coefficient  (Cd)  values  for  each  of  the  four  simulated  debris  elements  were  obtained 
by  measuring  the  flow  force  on  each  element  over  a range  of  velocities  in  a separate, 
fixed  boundary  flume.  The  fume  has  a width  of  1 .22  m,  and  a length  of  30.5  m. 
Approach  flow  velocities  were  measured  using  a pitot  tube  and  velocity  values  averaged 
over  a two  minute  sample  period.  Plate  4 shows  the  fixed  boundary  flume  and  force 
gauge  apparatus.  The  force  on  each  element  was  calculated  by  converting  the  average  of 
five  mass  readings  from  the  balance  to  a force  (N)  using  equation  46. 


Plate  4;  Force  gauge  set  up  in  fixed  boundary  flume  to  calibrate  drag  coefficient 
values  for  each  element. 

The  residual  value  caused  by  the  support  spindle  alone  was  subtraeted  from  each  reading 
and  the  actual  force  on  the  element  calculated  taking  a moment  about  the  spindle 
fulcrum.  Drag  force  on  eaeh  element  was  measured  over  a range  of  seven  depths.  Total 
flow  depth  during  the  test  runs  was  0.3  m.  The  results  are  presented  in  Table  8,  whieh 
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also  shows  minimum  and  maximum  balance  readings  at  eaeh  depth,  and  element 
Reynolds  and  Froude  numbers. 

Figure  9 shows  the  relationship  between  drag  eoefficient  and  flow  Reynolds  number 
(Re).  The  dark  curve  on  the  plot  shows  variation  of  Cd  with  Re  for  a smooth  infinitely 
long  cylinder  in  infinitely  wide  flow  (White,  1979).  For  bodies  with  a finite  length  drag 
coefficient  deereases  as  length  is  reduced  beeause  flow  round  the  ends  of  a body  reduees 
the  pressure  in  front  and  increases  the  pressure  behind  it.  Drag  coefficient  values  for 
eylinders  with  a range  of  1/d  ratios,  published  by  White  (1979),  are  shown  in  Figure  9 
(top  left  hand  box). 

Drag  results  for  the  four  elements  are  plotted  on  this  graph  as  eolored  symbols.  Ideally  all 
the  measured  results  should  cluster  about  the  1/d  = infinite  curve.  There  is  some  spread  of 
the  data  however,  both  above,  and  below  the  curve.  This  spread  can  be  explained  in  part 
beeause  1/d  ratios  vary  from  15  for  the  0.3m  element,  down  to  3.7  for  the  0.07m  element 
(see  box  in  bottom  left  corner  of  the  graph),  so  values  ean  be  expected  to  vary  between 
elements,  as  suggested  by  White's  (1979)  data.  White's  data  also  indieates  that  objects  of 
a finite  length  should  have  Cd  values  less  than  1.2.  However,  and  my  experiments 
produeed  Cd  values  as  high  as  1.5  (refer  to  Figure  9,  and  notiee  that  the  plotted  data 
extends  above  the  1/d  = infinite  curve).  We  attribute  this  diserepaney  to  the  fact  that  my 
experimental  flume  was  quite  narrow,  so  blockage  effects  (refer  to  section  2.2)  will  have 
been  signifieant,  whereas  White's  data  assumes  flow  of  infinite  width  (remember  that 
bloekage  increases  the  drag  eoefficient). 


45 


« 

.1  'i 


'VI  t I 


I 


Element  and 
velocity  depth 
from  bed  (m) 

Approach 

velocity 

(m/s) 

Balance 

reading 

(g) 

Min. 
reading 
(g) 

Max. 
reading 
(g) 

Residual 

reading 

Jg-) 

Drag 

force 

(N] 

Drag 

coeffieieni 

fn 

Reynolds 

number 

Re, 

Froude 

number 

Fr, 

0.3  m element  [element  area  (Ap)  = 0,0057m‘] 

029 

0219 

48  48 

47  32 

50,55 

19.18 

0.178322 

1.2801 

416  1 

0.507 

0.28 

0.218 

50  04 

48  7 

51  52 

19.59 

0 1 7633 

1 3018 

414,2 

0,504 

0,27 

0217 

46  82 

4627 

49.32 

20.01 

0 151 177 

1 1264 

412.3 

0 502 

0.26 

0.216 

43.94 

42,08 

46  98 

20.14 

0.130772 

0 9834 

410,4 

0 500 

0.25 

0 215 

42.08 

41.62 

45.09 

20,9 

0.1 13474 

08613 

408.5 

0 497 

0 2 

0.21 1 

42  97 

41.24 

47  16 

21 

0.108009 

0.8512 

400,9 

0,488 

0 15 

0.205 

44.83 

43.19 

47.93 

21.09 

0.101802 

0.8499 

389.5 

0.474 

0. 15  III  element  [element  area  (Ad)  = 0. 00285m’] 

0,29 

0,219 

27  1 1 

26,32 

27.72 

12,7 

0.085757 

1 254 

416  1 

0.507 

0.28 

0.218 

26.57 

25,66 

27  24 

14,62 

0.077988 

1 151 

414.2 

0.504 

0.21 

0.217 

27.37 

26.98 

28.01 

14,48 

0.072684 

1.083 

412,3 

0,502 

0.26 

0.216 

25.54 

25,19 

26.85 

13,88 

0.064067 

0.9636 

410.4 

0.500 

0.25 

0,215 

27  47 

27,7 

28,71 

15.57 

0,063755 

0.9678 

408.5 

04979 

0.2 

0.21  1 

27.95 

26.21 

27.3 

14.23 

0,054448 

0.8582 

400  9 

0.4887 

0.15 

0.205 

33  57 

32,7 

35  82 

17.26 

0,042239 

0,7053 

389.5 

04748 

0. 1 III  element  [element  area  (Ad)  = 0.001 9m'] 

0,29 

0.219 

21.07 

20,09 

21  78 

13.93 

0.069867 

1.5334 

416  1 

0.507 

0.28 

0 218 

22.34 

21,3 

22.59 

14.37 

0.061788 

1,3686 

4142 

0,504 

0.27 

0.217 

23.2 

22.6 

23.79 

15.39 

0.050915 

1,1382 

412.3 

0.502 

026 

0.216 

22,2 

21.95 

22.95 

14,21 

0 043902 

0.9904 

4104 

0.500 

0.25 

0,215 

24 

23.51 

24.19 

15.33 

0.04645 

1.0577 

408  5 

0497 

0.2 

0.21  1 

25  84 

24,88 

26 

16,35 

0 040681 

0 9618 

400  9 

0 488 

0 15 

0,205 

25.98 

24,87 

26.02 

16.62 

0.030403 

0 7613 

389.5 

0 474 

0.07m  element  [element  area  (Ad)  = 0.001 3m'] 

0,29 

0219 

18.83 

19  12 

19.98 

12  9 

0.035291 

1 1320 

416.1 

0,507 

0,28 

0.218 

19,81 

23  86 

20,38 

14,33 

0,031734 

1.0272 

414  2 

0.504 

0.27 

0,217 

18  5 

19.3 

19,36 

13.26 

0 029547 

0,9653 

412  3 

0,502 

0.26 

0,216 

18.02 

18.67 

22.7 

13,07 

0.027198 

0.8968 

410,4 

0 500 

0.25 

0,215 

21,7 

22,08 

23.34 

16  96 

0,025395 

0.8451 

408  5 

0.497 

0.2 

0,21 1 

22,72 

24.85 

23.25 

17 

0,023342 

0,8065 

400.9 

0 488 

0,15 

0.205 

21  35 

21,9 

22.02 

16.46 

0.020969 

0.7677 

389,5 

0,474 

where  force  (N)  = mass  (g)  / 1000  kg  x 9,81 , drag  coefficient  = force  / (p  x 0.5  x An  x v"),  element  Reynolds  number  Re^  = v d / v,  d 
element  diameter  (0,019  m),  and  element  Fronde  number  Fr^  = v / Vg  d 


Table  8:  Drag  coefficient  on  elements  at  a range  of  depths  in  infinite  flow 
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Figure  9:  Element  drag  coefficient  as  a function  of  Reynolds  number 


It  is  evident  from  Table  8 that  drag  coefficients  for  each  element  tend  to  increase  with 
distance  from  the  bed.  This  is  due  to  the  increase  in  flow  velocity  towards  the  flow 
surface.  This  relationship  becomes  more  evident  if  the  drag  coefficient  results  are  plotted 
as  a function  of  relative  flow  depth  (Figure  10).  Relative  flow  depth  is  defined  as  total 
flow  depth  (d)  divided  by  the  depth  from  the  center  of  the  element  to  the  channel  bed  (Y). 
Drag  coefficient  increases  slowly  with  relative  depth  below  Y/d  of  0.85.  However,  above 
this  level  drag  coefficient  increases  rapidly  to  a maximum  at  the  flow  surface.  This 
accelerated  increase  is  drag  close  to  the  water  surface  can  be  attributed  to  a phenomenon 
known  as  wave  drag.  This  occurs  when  a submerged  obstruction  in  flow  begins  to 
disturb  the  flow  surface,  setting  up  surface  waves.  Energy  is  required  to  create  these 
surface  waves  which  consequently  means  that  more  drag  is  being  imposed  by  the  object 
on  the  flow  than  the  simple  disturbance  of  stream  lines  around  the  object.  These  data 
show  that  drag  coefficients  may  be  significantly  underestimated  if  generally  estimated 
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coefficients  are  used  to  characterize  bluff  bodies  without  considering  the  objects 
proximity  to  the  fluid  free  surface. 


Figure  10:  Drag  coefficient  as  a function  of  relative  flow  depth 


7.2  Selection  of  drag  coefficient  values 

The  drag  coefficient  values  for  each  LWD  element  determined  in  the  fixed  boundary 
flume  were  then  used  to  solve  equation  47  to  determine  approach  flow  velocities  in  the 
Abiaca  Creek  model.  Values  selected  were  those  corresponding  to  the  appropriate 
element  depths  in  the  model  and  these  values  are  highlighted  in  bold  in  Table  8 and 
displayed  in  Table  9. 


Element  length  (m) 

Depth  of  flow  from 

Corresponding  drag  | 

surface  to  element 

coefficient  value  from 

center  (cm) 

Table  8 

0.07 

4.0 

0.84 

0.10 

4.0 

1.05 

0.15 

2.5 

1.08 

0.30 

0.5 

1.28 

Table  9:  Drag  coefficient  values  used  to  solve  equation  47 
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7.3  Estimation  of  friction  factor  due  to  LWD 

Shields  and  Gippel  (1995)  developed  a method  for  estimating  Darcy-Weisbach  friction 
factor  due  to  LWD  based  upon  the  work  of  Smith  and  Shields  (1992).  They  express  total 
friction  factor  as: 

fiotal  Igrain  IbedfonTi  & bar  Ibends  ^debris  (48) 

Friction  factor  due  to  grain  resistance  (fgram)  is  defined  by  equation  7 while  friction  factor 
due  to  bends  (fbends)  and  bedforms  and  bars  (fbedtorm  & fbar)  must  be  estimated.  The 
proportion  of  channel  roughness  ascribed  to  debris  is  detemiined  to  be: 


ZCo,A 


Di 


Cdebris  ^ 


(49) 


6ir  wL 

where,  a = kinematic  energy  correction  factor  (roughly  1.15).  and  i = debris  element  1,  2, 
etc. 


Equation  49  is  used  to  estimate  friction  factor  created  by  the  LWD  elements  in  the  model 
flume.  These  values  can  be  compared  with  measurements  reported  by  Shields  and 
Gippel  (1995)  from  data  taken  from  the  Obion  River  (Tennessee)  and  the  Tumut  River 
(Australia),  and  with  estimates  of  friction  factor  for  the  prototype  reach  on  Abiaca  Creek 
reported  by  Watson  et  al.  (1993). 

Drag  coefficients  used  to  solve  equation  47  were  those  presented  in  Table  9.  Results  are 
shown  in  Table  10.  Comparison  of  these  results  with  those  presented  by  Shields  and 
Gippel  (1995)  (see  Table  1 1)  demonstrates  that  values  calculated  for  the  flume  model  are 
of  the  same  order  of  magnitude  as  those  estimated  in  the  field  situation. 

The  flume  model  can  therefore  be  regarded  to  have  accurately  re-created  the  average  drag 
and  friction  factor  properties  of  LWD  as  found  in  the  field. 

Total  friction  factor  in  the  prototype  reach  on  Abiaca  Creek  is  suggested  by  Watson  et  al. 
(1993)  to  be  n = 0.032.  Converting  this  value  to  a Darcy-Weisbach  friction  factor  using 
the  following  equation  (Richards,  1982): 
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f=8gn“/R'^^  = 0.235 


(50) 


Element 

Length 

(m) 

Cd  in 
infinite 
flow 

Element  area 
normal  to  flow 
(m') 

fdebris 

Ftotai 

0.30 

1.28 

0.0057 

0.051 

0.146 

0.200 

0.15 

1.08 

0.0029 

0.022 

0.146 

0.168 

0.10 

1.05 

0.0019 

0.014 

0.146 

0.160 

0.07 

0.84 

0.0013 

0.008 

0.146 

0.154 

Table  10:  Calculation  of  Darcy-Weisbach  friction  factor  due  to  LWD 


Site 

Q 

(mVs) 

V 

(m/s) 

s 

(ni/km) 

R 

(m) 

D51, 

(mm) 

ff>cd 

Abends 

^debris 

f,„u, 

Obion  River 

3 9 

0,29 

0 63 

0.9 

0 27 

0.22 

0,00 

0.20 

0,42 

Obion  River 

4.0 

0.33 

0.61 

0.7 

0 44 

0,16 

0,00 

0.18 

0.34 

Obion  River 

90 

0.46 

0 78 

II 

0.44 

0,17 

0.00 

0.17 

0.34 

Obion  River 

16,2 

0,53 

0.45 

16 

0 27 

0 1 1 

0,00 

0.04 

0.16 

Obion  River 

16.9 

0.65 

0,51 

14 

0,44 

0.09 

0,00 

0,02 

0,1  1 

Obion  River 

36.9 

0.72 

0.46 

2.7 

0.27 

Oil 

0.00 

0,07 

0,18 

Obion  River 

40  5 

0 86 

0 60 

2.5 

0 44 

0.10 

0 00 

0.03 

0,13 

Obion  River 

409 

0.77 

0 49 

2.7 

0,27 

0,10 

0,00 

0.07 

0 17 

Obion  River 

41  3 

0.86 

0 59 

2.5 

0 44 

0,10 

0 00 

0,03 

0.13 

Obion  River 

3 6 

0 46 

0.42 

0.5 

0 44 

0.06 

0.00 

0.03 

0.1  1 

Obion  River 

4 2 

0 45 

0.54 

0,6 

0,57 

0,09 

0.00 

0.00 

0.09 

Obion  River 

10.7 

0.57 

0 42 

1.2 

0.57 

0.08 

0.00 

0,06 

0.14 

Obion  River 

22.2 

0 84 

0.46 

16 

0.57 

0.06 

0 00 

0,01 

0.07 

Obion  River 

47,9 

0,91 

0 52 

2.7 

0.57 

0.09 

0.00 

0,01 

0 10 

Obion  River 

53.2 

098 

0.54 

2 8 

0 57 

0.08 

0.00 

0.01 

0 10 

Tumut  River 

112 

1,01 

1 10 

2.8 

16.0 

0.1 1 

0,04 

0,02 

0,17 

Tumut  River 

96 

0 90 

1.00 

2.9 

22.0 

0.13 

0.07 

0,03 

0.22 

Tumut  River 

96 

0 90 

1,50 

29 

28.0 

0,14 

0.12 

0.00 

0.26 

Tumut  River 

125 

1.08 

1,30 

3,0 

16,0 

0,10 

0,04 

0,00 

0.14 

Tumut  River 

92 

0.93 

0.90 

2,6 

22,0 

0.13 

0.06 

0,00 

0 19 

Tumut  River 

92 

0.93 

1,50 

2.6 

28.0 

0,14 

0,1 1 

0.00 

0 25 

Table  11:  Computation  of  Darcy-Weisbach  friction  factor  in  the  Obion  River  and 


Tumut  River  (after  Shields  and  Gippel,  1995) 


This  prototype  value  lies  outside  the  upper  range  of  ftotai  values  predieted  in  Table  10. 
Considering  that  the  prototype  reach  contains  bed  features  that  are  incorporated  into 
friction  factor  in  equation  50  but  were  not  present  in  the  model  situation,  it  would  appear 
that  Watson’s  field  estimate  of  total  friction  factor  is  somewhat  low.  This  may  be  due  to 
neglecting  the  important  influence  of  LWD. 

From  the  model  flume  results,  and  those  presented  by  Shields  and  Gippel  (1995)  it  may 
be  concluded  that  friction  factor  in  natural  channels  may  well  be  considerably 
underestimated.  This  is  because  even  if  resistance  caused  by  bends  and  bedforms  is 
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taken  into  account,  reach  averaged  resistance  due  to  debris  accumulations  is  rarely 
considered  but  can  create  as  much  frictional  resistance,  if  not  greater  than  the  resistance 
offered  by  channel  bed  roughness  (see  Table  1 1 ). 

7.4  Drag  force  on  LWD  elements  in  the  model 

Drag  force  on  each  element  is  plotted  against  time  in  Figure  1 1.  The  temporal  change  in 
force  readings  can  be  explained  by  the  fact  that  as  an  element  was  lowered  into  the  model 
channel,  flow  was  forced  to  accelerate  past  it  in  order  to  maintain  continuity  of  discharge. 
The  affect  of  flow  acceleration  and  increased  turbulence  was  to  raise  boundary  shear 
stress  above  the  threshold  condition  causing  the  channel  bed  and  banks  to  scour  in  the 
vicinity  of  the  element.  The  scoured  sediment  was  carried  downstream  and  deposited 
where  flow  decelerated  to  the  point  of  losing  competence  to  maintain  transport. 
Flowever,  the  rate  of  scouring  declined  over  time  because  as  the  cross-sectional  area  of 
the  channel  increased  so  velocities  fell  again  thus  reducing  boundary  shear  stress.  This 
negative  feedback  between  channel  enlargement  and  boundary  shear  stress  reduction 
eventually  led  to  an  equilibrium  state  once  more,  where  the  channel  boundary  returned  to 
a threshold  state,  but  with  a new,  enlarged  cross-sectional  area.  Flow  velocity  around  the 
element  therefore  decayed  in  an  asymptotic  manner  through  time,  causing  the  drag  force 
experienced  by  the  element  to  decline  in  a similar  manner. 

Power  functions  are  often  found  to  best  characterize  changes  in  energy  and  landscape 
adjustment  in  geomorphic  systems  (Thorn,  1988).  This  fact  gave  the  justification  for 
fitting  power  functions  to  the  model  data,  and  regression  analysis  demonstrates  that 
power  decay  functions  describe  the  measured  changes  in  drag  force  very  well.  The 
equations  of  these  functions  and  their  respective  r‘  values  are  shown  in  Figure  1 1 and  in 
Table  12. 
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Figure  11:  Measured  change  in  drag  force  over  time  on  each  LWD  element  in  the 
flume  model 


Element  Length 
(m) 

Equation 
constant  (a) 

Equation 
Exponent  (b) 

Goodness  of  fit 
(r') 

0.30 

0.213 

-0.048 

0.94 

0.15 

0.132 

-0.034 

0.95 

0.10 

0.044 

-0.075 

0.92 

0.07 

0.004 

-0.249 

0.97 

Table  12:  Equations  describing  drag  force  decay  on  LWD  elements.  Regression 
Functions  take  the  form:  Force  = a + Time'’ 


It  is  evident  from  Figure  1 1 and  Table  12  that  the  larger  the  element  the  greater  the  total 
force  acting  upon  it. 


7.5  Drag  force  on  LWD  elements  in  the  prototype 

Figure  12  shows  the  change  in  drag  force  over  time  on  LWD  elements  in  the  prototype 
reach  (Abiaca  Creek)  computed  from  the  model  data.  The  distribution  of  the  data  is 
obviously  the  same  as  for  the  model,  but  absolute  values  have  been  expanded  according 
to  the  force  and  time  scaling  relationships  described  in  section  7.  Note  that  drag  force  is 
measured  in  kilonewtons  and  that  the  time  scale  extends  up  to  6 hours.  Prototype  debris 
element  dimensions  shown  in  Figure  12  are  calculated  in  Table  7. 
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Power  functions  have  been  fitted  to  the  data  and  both  the  equations  and  values  are 
shown  in  Table  13.  The  shape  of  each  fitted  functions  is  obviously  the  same  as  for  the 
model  data,  the  actual  values  are  simply  larger  in  accordance  with  the  scaling  ratios. 


Prototype  Iilement 
Length 

4 18m  (Underflow  Jam) 
9m  (Dam  Jam) 

A 6m  (Deflector  Jam) 

• 4.2m  (Parallel  Jam) 


Figure  12:  Change  in  drag  force  over  time  on  LWD  elements  in  the  prototype  reach 
computed  from  model  data 


Element  Length 
(m)/  Jam  type 

Equation 
constant  (a) 

Equation 
Exponent (b) 

Goodness  of  fit 
(••') 

18m 

21.059 

-0.048 

0.94 

9m 

12.987 

-0.034 

0.95 

6m 

4.400 

-0.075 

0.93 

4.2m 

0.515 

-0.241 

0.98 

Table  13:  Equations  describing  drag  force  decay  on  LWD  elements  in  the  prototype. 
Regression  Functions  take  the  form:  Force  = a + Time'* 


For  example,  we  can  expect  in  the  prototype  river  during  a Q2  discharge  event  that  a tree 
trunk  with  a length  of  18  m and  a diameter  of  1.1  m will  experience  a drag  force  of  25.5 
kN  when  it  falls  into  the  channel,  normal  to  the  flow,  and  that  this  drag  force  will 
decrease  over  a period  of  6 hours  to  17.5  kN.  If  one  assumes  that  total  force  decay  time 
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is  equivalent  to  the  time  it  takes  the  channel  morphology  to  adjust,  the  geomorphic 
changes  to  a sand-bed  creek  caused  by  LWD  elements  are  complete  after  a short  period 
of  time.  Obviously,  if  a tree  were  to  fall  into  such  a creek  while  flow  conditions  were 
less  than  the  2-year  event  total  force,  adjustment  times  would  be  longer  by  an  amount 
proportional  to  the  approach  flow  velocity  for  that  discharge. 

It  can  be  concluded  that  the  range  and  temporal  change  in  values  in  Abiaca  Creek  caused 
by  the  channel  boundary  adjustment  would  be  very  similar  to  those  predicted  by  the 
flume  model.  Given  that  the  prototype  reach  was  chosen  as  representative,  these  results 
must  also  be  more  generally  applicable  to  LWD  interaction  with  the  flow  in  any  sand-bed 
river  in  Northern  Mississippi. 

7.6  Incident  flow  velocity  on  LWD  elements  in  the  model 

Change  in  incident  velocity  on  each  element  over  time  was  computed  by  substituting  the 
measured  drag  force  data  (Figure  1 1 ) and  drag  coefficient  values  (Table  8)  into  equation 
47.  The  results  are  displayed  in  Figure  13.  Given  that  all  other  variables  in  equation  47 
are  constant  for  each  time  step  except  drag  force,  it  is  not  surprising  to  find  that  incident 
flow  velocity  on  each  element  decays  over  time  in  exactly  that  same  manner  as  for  the 
drag  force  measurements.  The  data  were  fitted  by  power  decay  functions,  and  function 
constants,  exponents,  and  goodness  of  fit  (r  values)  for  each  element  are  displayed  in 
Table  14. 
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Figure  13:  Change  in  approach  velocity  over  time  on  LWD  elements  in  the  flume 
model 


Regression  Functions  take  the  form:  Velocity  = a + Time*^ 

Element  Length 

Equation 

Equation 

Goodness  of  fit 

(m) 

constant  (a) 

Exponent (b) 

(•*') 

0.3 

0.342 

-0.024 

0.94 

0.15 

0.296 

-0.017 

0.95 

0.1 

0.210 

-0.037 

0.92 

0.07 

0.094 

-0.124 

0.97 

Table  14:  Equations  describing  incident  velocity  decay  at  each  LWD  element  in  the 
model 


It  might  be  expected  that  these  velocity  curves  would  converge  towards  an  asymptotic 
value  over  time  as  an  equilibrium  state  is  reached  in  the  channel  through  the  course  of 
channel  adjustment.  It  is  also  surprising  that  two  curves  in  Figure  13  fall  below  the  mean 
flume  flow  velocity  of  0.24  m/s  because  velocities  near  the  elements  must  have  been 
considerably  higher  in  order  to  cause  sediment  transport.  These  anomalies  can  be 
explained  by  the  fact  that  elements  were  not  suspended  at  the  same  depth  in  the  model. 
The  0.3  m element  was  held  Just  below  the  flow  surface  while  the  0.15  m,  0.1  m and  0.07 
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m elements  were  suspended  at  progressively  greater  depths  in  the  flow.  Larger  elements 
will  have  experienced  higher  approach  flow  velocities  than  smaller  elements,  regardless 
of  the  disturbance  caused  by  the  element  itself. 

Figure  13  clearly  demonstrates  that  the  increase  in  flow  velocity  caused  by  the  elements 
reduction  of  the  channel  cross  sectional  area  is  quickly  accommodated  by  the  channel. 
This  occurred  through  erosion  of  the  channel  boundary  causing  velocity  to  fall  in  an 
asymptotic  manner  over  time  until  boundary  shear  stress  returned  to  the  critical  value  and 
sediment  transport  ceased. 


7.7  Incident  flow  velocity  on  LWD  elements  computed  for  the  prototype 

Figure  14  shows  the  change  in  approach  flow  velocity  for  the  four  different  jam  types  as 
predicted  for  the  prototype.  Power  decay  functions  have  also  been  fitted  to  these  data. 
The  fitted  equation  constants,  exponents  and  goodness  of  fit  r^  values  are  shown  in  Table 
15. 


Prototype  Element 
Length 

♦ 18m  (Underflow  Jam) 

• 4.2m  (Parallel  Jam) 

A 6m  (Deflector  Jam) 

9m  (Dam  Jam) 


Figure  14:  Change  in  LWD  element  approach  flow  velocity  over  time  in  prototype 
reach  computed  from  model  data 
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Element  Length 
(m)  / Jam  Type 

Equation 
Constant  (a) 

Equation 
Exponent  (b) 

Goodness  of 
Function  Fit 

(r^) 

18m 

1.801 

-0.024 

0.94 

9m 

1.553 

-0.017 

0.95 

6 m 

1.113 

-0.038 

0.93 

4.2m 

0.515 

-0.125 

0.98 

Table  15:  Equations  describing  incident  velocity  decay  on  LWD  elements  in  the 
prototype.  Regression  Fnnctions  take  the  form:  Velocity  = a + Time^ 


The  behavior  of  velocity  change  for  each  jam  type  is  the  same  as  that  discussed  for  the 
model  data  but  actual  values  are  larger.  Thus,  maximum  total  adjustment  time  in  the 
prototype  is  6 hours.  One  can  expect  during  a Q2  discharge  event  that  a tree  trunk  with  a 
length  of  18  m and  a diameter  of  1 .1  m will  experience  an  incident  flow  velocity  of  0.385 
m/s  when  it  falls  into  the  channel,  normal  to  the  flow,  and  that  this  velocity  will  decrease 
over  a period  of  6 hours  to  0.338  m/s. 


7.8  Conclusions  from  the  hydraulic  data 

The  prototype  drag  force  and  incident  velocity  results  lead  to  the  conclusion  that  LWD 
obstructions  are  quickly  accommodated  for  by  geomorphic  adjustment  in  mobile  sand- 
bed  channels.  So,  even  dense  LWD  jams  in  sand-bed  rivers  will  not  disturb  the  reach- 
scale  flow  regime  to  the  extent  that  general  scour  (i.e.  bed  degradation)  might  be 
triggered.  Only  local  bed  and  bank  adjustments  are  likely.  It  can  be  further  concluded 
that  LWD  in  sand-bed  channels  that  are  undergoing  degradation  and  lateral  adjustment  as 
described  by  the  Schumm  et  al.  (1984)  Channel  Evolution  Model  is  unlikely  to  act  as  a 
significant  structural  barrier  to  stop  or  slow  the  upstream  migration  of  a knick-point.  In 
short,  the  flow  will  very  quickly  erode  both  the  bed  and  banks  in  the  immediate  vicinity 
of  a jam  to  bypass  the  obstruction.  LWD  in  sand-bed  rivers  is  unlikely  to  alter  the  gross 
long-term  energy  balance  (the  ratio  of  potential  to  kinetic  energy)  of  a reach  by 
significantly  damming  up  the  flow,  a finding  which  is  a direct  contrast  to  that  reported  by 
Marston  (1982)  for  LWD  jams  in  gravel  bed  rivers. 

The  geomorphic  features  associated  with  LWD  jams  in  unstable  sand-bed  rivers  are 
assumed  to  become  relatively  stable  very  quickly  after  the  debris  has  entered  the  channel. 


57 


1 


rit  .1 


I 


I 


I 


It 


1 


» , 


^ I 


These  conclusions  agree  with  field  observations  as  the  local  channel  geomorphology, 
such  as  upstream  and  downstream  bars  and  the  depth  and  extent  of  scour  holes,  at  the 
majority  LWD  Jams  monitored  over  the  two  year  survey  period  varied  very  little  over 
time.  Indeed,  local  bank  failure  scars  and  in-channel  bars  caused  by  the  six  new  jams 
that  had  formed  in  the  survey  reaches  between  the  1995  and  1996  were  found  to  have 
already  been  colonized  by  a dense  vegetation  of  grasses  and  shrubs  by  the  time  of  the 
second  survey. 

7.9  Geomorphic  results:  Flume 

The  previous  discussion  has  focused  on  the  hydraulic  effects  of  the  LWD  elements.  The 
next  step  is  to  examine  how  these  effects  interact  which  the  channel  boundary  to  produce 
a set  of  LWD  induced  geomorphic  features. 

At  the  end  of  each  test  run,  once  all  sediment  transport  had  been  observed  to  cease,  the 
discharge  in  the  tlume  was  gradually  reduced  until  all  flow  drained  into  the  sand 
underneath  the  cut  channel.  The  flume  was  then  surveyed  along  a length  from  0.1  m 
upstream  of  the  LWD  element  position  to  a point  where  channel  topography  was  no 
longer  disturbed  downstream  of  the  elements.  A point  gauge,  mounted  on  the  instrument 
earriage,  was  used  to  measure  break  of  slope  across  the  ehannel  at  regular  distances  down 
the  flume,  which  thus  produced  a set  of  x,y,z  data  describing  the  channel  topography. 
This  data  was  then  entered  into  a surfaee  interpolation  program  that  generated  eontour 
maps  of  the  channel.  The  channel  was  also  surveyed  prior  to  each  test  run  and  a surface 
generated.  These  data  enabled  the  generation  of  contour  maps  of  relative  channel 
adjustment  by  subtraeting  the  pre-test  surface  from  the  post-test  surface.  The  resulting 
contour  surfaces  for  each  test  are  shown  in  Figures  15  and  16.  The  uneroded  channel 
‘floodplain’  between  the  channel  bank  top  and  the  flume  sidewalls  is  colored  in  gray 
(marked  left  and  right  banks),  without  a contour  surface  in  order  to  avoid  confusion  over 
the  channel  limits. 

In  Figures  15  and  16,  there  are  the  two  bands  of  aggradation  running  in  parallel  strips 
along  the  channel  at  about  0.2  and  0.45  m across  the  channel.  These  bands  with  heights 
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of  roughly  +40  mm  represent  upper-bank  material  that  has  been  dislodged  by  surface 
flow  disturbances  slumped  to  the  bank  base.  The  upper  bank  regions  show  this  sediment 
removal  and  slumping. 

A second  set  of  features  of  great  significance  are  the  waves  of  upper  bank  scour,  seen  on 
the  floodplain  edge  that  become  more  pronounced  the  larger  the  element.  These  waves 
of  scour  are  not  related  to  the  zones  of  bed  and  bank  scour  caused  by  flow  impingement 
on  the  channel  bed  and  lower  bank.  Instead  they  are  caused  by  flow  surface  waves 
hitting  the  channel  banks.  These  waves  are  generated  by  the  combined  effect  of  surface 
gravity  waves  set  up  by  flow  acceleration  under  the  elements  and  by  the  closure  of 
pressure  generated  streamlines  around  the  element  ends  (White,  1994).  The  form  and 
effect  of  the  waves  are  shown  schematically  in  Figure  17.  The  closer  an  object  is  to  the 
flow  surface,  the  more  pronounced  the  water  waves.  This  explains  why  the  0.15  m and 
0.3  m elements  (which  lay  at  0.20  and  0.05  m below  the  free  surface)  have  created  more 
pronounced  nodes  of  scour  than  the  0.07  m and  0.1  m elements  (which  both  lay  at  0.35  m 
below  the  free  surface).  Plate  5 shows  the  topographic  results  of  the  0.15  m element. 
The  first  set  of  pressure-generated  waves  on  the  downstream  side  of  the  element  are 
elearly  visible,  as  is  the  zone  of  element  induced  flow  disturbance,  that  diminishes  with 
distance  downstream.  This  contrasts  with  the  calm  flow  conditions  approaching  the 
element  (left-hand  side  of  the  picture).  The  largest  volumes  of  bank  erosion  were  caused 
by  deflection  of  subsurface  flow  due  to  flow  acceleration  around  the  LWD  elements. 
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Figure  15:  Contour  plots  of  relative  depth  change  for  the  flume  model:  0.07  m and  0.1  m debris  elements 
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Figure  16:  Contour  plots  of  relative  depth  change  for  the  flume  model:  0.15  m and  0.3  m debris  elements 
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Plate  5:  Flow  disturbance  created  by  the  0,15  m debris  element 
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It  is  worth  noting  that  the  bank  tops  were  located  at  0. 1 5 m and  0.53  m across  the  channel 
from  the  left  flume  side-wall  prior  to  element  immersion.  All  four  elements  can  therefore 
be  seen  to  have  caused  floodplain  failure  along  the  length  of  the  survey  reach. 

The  0.1  m (deflector)  debris  element  created  the  greatest  local  bank  erosion  (almost  0.1 
m into  the  right  floodplain)  because  it  diverted  a large  part  of  the  central  fast  flowing 
water  from  left  to  right  across  the  channel.  The  resulting  bank  failure  contrasts  greatly 
with  the  limited  bank  erosion  generated  by  the  0.07  m (flow  parallel)  element,  which  is 
only  0.03  m shorter,  but  lay  closer  to  the  base  of  the  left  bank  and  therefore  did  not 
disturb  the  central  fast  flowing  stream.  The  0. 1 5 m (Dam)  element  split  the  flow  into  two 
streams  that  impinged  on  each  bank  downstream  of  the  element.  This  caused  scour  of  the 
left  bank  to  a maximum  depth  of  12  mm  net  erosion,  and  a more  elongated  but  narrower 
band  of  erosion  along  the  right  bank  to  a maximum  depth  of  15  mm  net  erosion.  The 
eroded  material  was  transported  downstream  across  the  bank  face  and  deposited,  as 
lateral  bars  in  the  zone  of  low  shear  stress  that  occurs  at  the  break  of  slope  between  the 
bank  and  bed.  The  0.3  m element  lay  close  to  the  flow  surface,  spanning  the  full  width  of 
the  channel.  This  meant  that  flow  was  deflected  underneath  the  element  but  also  around 
either  end  which  produced  two  zones  of  high  shear  stress,  eroding  the  banks  up  to  a depth 
of  0.02  m of  net  change. 

The  location  and  severity  of  erosion  depends  on  the  size  of  the  debris  element  but  more 
importantly  upon  its  position,  both  vertically  in  the  flow  stream  and  laterally  across  the 
channel,  relative  to  the  core  of  maximum  flow  velocity  (White,  1994).  The  degree  of 
flow  impingement  on  either  bank  was  greatest  for  the  0.15  m and  0.1  m elements  for  this 
very  reason.  In  contrast  the  0.07  m element  created  little  bank  failure  as  it  did  not 
interfere  with  the  main  flow  stream.  While  the  0.3  m element  caused  the  greatest  net 
bank  erosion  in  small  zones  at  either  end  of  the  element,  erosion  downstream  was 
limited.  This  is  flow  was  diverted  mainly  under  the  obstruction  parallel  to  the  banks, 
rather  than  laterally  across  the  channel  that  set  up  a back  and  forth  oscillation  of  the  flow 
field  in  the  downstream  reach. 
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I'he  second  component  of  geomorphic  change  was  adjustment  of  the  channel  bed 
topography  by  the  disturbed  flow  field  (refer  to  Figures  15  and  16).  The  0.07  m element 
created  a scour  hole  in  the  center  of  the  channel  that  was  8 mm  deep  at  its  maximum. 
The  sand  scoured  from  this  hole  was  deposited  mainly  in  the  zone  of  calm  water  in  the 
wake  of  the  element,  forming  a lateral  bar  14  mm  high  and  0.3  m long  against  the  left 
bank.  A second  bar  formed  at  the  base  of  the  opposite  bank,  slightly  further  downstream, 
where  the  flow  turbulence  diminished.  Flow  then  oscillated  very  slightly  between  these 
two  features.  This  pattern  of  sedimentation  resembles  a set  of ‘alternate  bars'  (Richards, 
1982)  that  are  common  features  in  straight  sections  of  both  sand  and  gravel  bed  rivers.  A 
second  zone  of  bed  scour  developed  downstream  of  the  alternate  bars.  This  secondary 
scour  feature  was  observed  to  occur  downstream  of  the  zone  of  sediment  deposition  and 
is  thought  to  have  been  caused  by  flow  acceleration  over  the  crest  of  and  between  the  bar 
formations. 

The  0.1  m element  produced  a large  bar  in  its  wake,  which  was  0.35  m long  and  13  mm 
high.  The  material  that  built  this  bar  came  mainly  from  the  area  of  intense  bank  erosion, 
slightly  upstream  on  the  opposite  bank,  and  was  carried  across  the  channel  by  high- 
energy,  sediment-laden  flow  to  the  area  of  low  energy  in  the  element's  lee.  A second, 
elongated  but  lower  bar  formed  against  the  base  of  the  right  bank  running  between  0.05 
m and  0.75  m along  the  survey  reach.  A third  shorter  bar  formed  against  the  base  of  the 
left  bank  at  the  very  downstream  end  of  the  survey  reach.  This  pattern  of  alternating  bars 
demonstrates  that  the  element  set  up  an  oscillatory  flow  pattern  in  the  downstream  reach. 
The  erosive  strength  of  this  meandering  flow  was,  however,  dissipated  after  the  initial 
deflection  to  the  right  bank.  It  is  therefore  evident  that  such  an  obstruction  could  not 
disturb  the  flow  field  to  the  extent  that  reach  scale  channel  meandering  might  be  induced. 
Bed  scour  also  occurred  underneath  the  element  but  the  extent  and  depth  of  this  feature 
was  not  as  great  as  that  created  by  the  0.07  m element. 

The  0.15  m element’s  position  in  the  center  of  the  channel  caused  a large  body  of  the 
flow  to  be  diverted  downwards  onto  the  channel  bed  producing  a large  scour  hole  0.4  m 
long  and  up  to  8 mm  deep.  The  sediment  scoured  from  this  hole  was  deposited  as  a low 
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mid-channel  bar.  Flow  acceleration  over  this  feature  generated  a second  zone  of  bed 
scour  at  the  lower  end  of  the  survey  reach.  This  set  of  bedforms  closely  resembled  a 
elassic  pool-riffle  sequence  (Knighton,  1984),  although  there  was  no  obvious  sediment 
grain  size  sorting  in  the  model. 

The  0.3  m element  created  the  most  extensive  and  deepest  scour  hole  (12  mm  of  net 
scour)  because  it  blocked  and  subducted  the  fastest  moving  body  of  flow  near  the  free 
surfaee.  The  scoured  sediment  was  then  deposited  as  a large  mid-channel  bar,  which 
extended  across  the  full  width  of  the  channel  bed  up  to  a maximum  height  of  14  mm. 
Flow  aceeleration  over  this  bar  produeed  an  equally  large  zone  of  secondary  bed  seour 
towards  the  downstream  end  of  the  survey  reach.  The  channel  adjustments  eaused  by  the 
0.3  m element  are  shown  in  Plate  6. 


Plate  6:  Channel  adjustment  caused  by  the  0.3  m debris  element 

The  extent  and  depth  of  bed  scour  is  therefore  seen  to  be  a function  of  debris  element  size 
and  also  position  with  respeet  to  the  velocity  profile  (both  vertieally  from  the  bed  and 
laterally  across  the  ehannel).  Large  elements  and/or  elements  near  to  the  center  of  the 
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channel  and  high  in  the  velocity  profile  cause  a greater  depth  and  extent  of  scour  than 
smaller  elements  and/or  elements  near  to  the  channel  bed  and  banks.  The  size  of  mid 
channel,  alternate  and  lateral  bars  is  proportional  to  the  volume  of  material  eroded  from 
the  scour  holes. 

7.10  Geomorphic  results:  Prototype 

Because  the  flume  is  a fully  scaled  model  the  geomorphic  results  can  be  scaled  to 
prototype  dimensions,  and  compared  directly  with  the  field  results,  rather  than  having  to 
resort  to  comparison  by  ‘similarity  of  process'  (Cherry  and  Beschta,  1986;  Young,  1991; 
Peakall  et  ak,  1996;  Braudrick  et  ak,  1997).  The  fume  survey  data  was  therefore  scaled 
to  prototype  dimensions  using  the  X and  Y scaling  relationships  defined  by  the  model  (Xr 
= 59.667;  Yr  = 27.142).  The  rescaled  data  were  then  converted  to  contour  plots.  The 
resulting  channel  topographies  are  shown  in  Figures  18  and  19.  Note  the  difference  in 
scale  between  Figures  15  and  16  (lengths  and  depths  in  cm)  and  Figures  18  and  19 
(lengths  and  depths  in  m).  The  reach  length  surveyed  from  the  debris  element 
downstream  is  60m,  and  the  undisturbed  pre-test  channel  top  width  is  17.9  m.  These 
surfaces  represent  the  actual  scale  of  features  that  would  be  found  in  the  prototype  reach 
of  Abiaca  Creek  caused  by  each  debris  element  type  after  6 hours  of  immersion  in  the 
flow  during  the  2-year  discharge  event.  However,  the  key  premise  behind  this 
experimentation  is  that  the  geomorphology  associated  with  LWD  jams  varies  in  a 
predictable  manner  downstream  though  the  channel  network  as  channel  width  increases. 
Channel  width  increases  in  a predictable  manner  (through  the  use  of  regime  equations)  or 
can  be  estimated  through  the  knowledge  of  other  important  processes  (channel  evolution 
models)  with  distance  downstream  through  the  drainage  basin. 

As  long  as  the  average  ratio  between  channel  width  and  debris  length  that  characterizes 
each  jam  type  was  reproduced  by  varying  debris  length,  the  suite  of  geomorphic  forms 
generated  should  be  the  same  as  those  found  in  the  field.  Although  the  model  had  to  be 
scaled  from  a particular  field  site,  the  hydraulic  and  geomorphic  results  generated  can  be 
regarded  to  be  applicable  to  any  river  that  has  a mobile  sand  bed  and  highly  erodible 
banks  and  similar  morphometric  and  hydraulic  characteristics. 
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The  surfaces  displayed  in  Figures  18  and  19  have  exactly  the  same  set  of  features  as 
those  in  Figures  15  and  16  so  it  is  not  necessary  to  describe  their  characteristics  and  mode 
of  formation.  The  prototype  surfaces  have,  however,  been  redrawn  in  a simplified  form 
in  Figures  20  and  21,  to  show  the  spatial  and  vertical  dimensions  of  bars  and  scour  holes 
that  can  be  attributed  solely  to  simple  LWD  forms  in  the  field. 
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Flow  direction  is  from  right  to  left 

Figure  18:  Contour  plots  of  relative  depth  change  calculated  for  the  prototype  reach:  Flow  Parallel  and  Deflector  Debris 
elements 
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Figure  19:  Contour  plots  of  relative  depth  ehange  calculated  for  the  prototype  reach:  Dam  and  Underflow  debris  elements 
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Figure  20:  Geomorphic  adjustments  to  a simple  trapezoidal  ehannel  caused  by  4.2  and  6-m  long  LWD  elements. 
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8 Conclusions 

The  geomorphic  results  displayed  in  Figures  20  and  21  answer  a number  of  questions  that 
have  arisen  through  the  course  of  this  research.  The  results  show  that  LWD  has  a distinct 
geomorphic  impact  upon  highly  erodible  channels,  causing  both  bed  scour  and  bank 
erosion,  and  the  formation  of  mid-channel,  lateral  and  alternate  bars. 

Secondly,  the  flume  results  also  help  to  solve  the  ‘chicken-and-egg'  problem  of  cause 
and  effect  which  arises  from  simple  observational  geomorphology.  In  this  case,  the 
problem  of  whether  debris  jams  cause  pool-riffle  formations  or  whether  debris  simply 
collects  where  it  runs  aground  against  the  upstream  face  of  riffles  and  bars.  The 
geomorphic  results  clearly  demonstrate  that,  whilst,  in  some  instances  debris  may  be 
trapped  by  high  points  in  the  channel  bed  topography,  it  will  also  create  a distinctive  set 
of  geomorphic  features  regardless  of  pre-existing  channel  forms. 

Given  that  the  test  runs  continued  until  sediment  erosion  and  transport  ceased,  it  is 
interesting  to  find  that  the  impact  of  debris  elements  is  local  in  nature  as  the  channel 
becomes  progressively  less  disturbed  with  distance  downstream  from  each  element.  Of 
the  three  important  types  of  scour  (Local,  Constriction,  and  General;  see  Raudkivi; 
1990),  only  local  scour  and  constriction  scour  are  likely  at  jams.  General  scour  and 
degradation  were  not  observed  in  the  flume  tests.  With  regard  to  the  question  of  channel 
planform  adjustment  caused  by  debris,  it  is  evident  from  the  flume  data  that  even  the 
deflector  element  did  not  to  cause  lateral  oscillation  of  the  flow  to  the  e.xtent  that 
meanders  developed  in  the  downstream  reach. 

These  conclusions  are  based  solely  on  the  flume  results  and  the  limitations  of  this  model 
must  be  considered.  There  was,  for  example,  no  sediment  transport  from  upstream,  so 
backwater  sediment  wedges,  which  were  observed  to  form  upstream  of  dam  type  jams  in 
the  field,  could  not  form  in  the  model.  The  model  ‘jams'  were  composed  of  single 
elements  whereas  ‘real’  jams  are  usually  more  complex  features  so  the  extent  of  the 
geomorphic  adjustment  in  the  model  may  be  considerably  less  than  the  adjustment  that 
might  be  cause  by  large  and  compact  (watertight)  jams  in  the  field.  Indeed,  some  authors 
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have  reported  debris  jams  that  have  caused  channel  avulsion  (Keller  and  Swanson,  1979). 
Because  the  model  bank  materials  were  composed  only  of  sand,  it  is  not  known  whether 
the  extent  of  bank  erosion  would  be  greater  (further  into  the  floodplain)  or  less  had  bank 
failure  taken  place  through  both  fluvial  erosion  and  by  mass  failure.  However  the  results 
presented  herein  represent  the  simplest  geomorphic  forms  that  the  field  scientist  can 
expect  to  see  in  association  with  the  different  jam  types  encapsulated  in  the  Debris  Jam 
Classification  Model. 
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APPENDIX 


Large  Woody  Debris  Formation  Survey  used  by  Smith  and  Shields  (1992) 


Stream  Name 

Reach 

Information 


Date 


.Time. 


Width-Perpendicular  to  Flow  Direction 


L<B/2 


B/2<L<B 


W<B/4 


B/4<W<B'/2 


~B/2^W^ 


L>B 

Length- 
Parallel 
to  Flow 
Direction 


TYPE  A:  COLLAPSED  BRIDGE 


L<B/2 


B/2<L<B 


L>B 


Length- 
Parallel 
to  Flow 
Direction 


TYPE  C:  DRIFT 


Width-Perpendicular  to  Flow  Direction 


W<B/4 


B/4<W<b/2 


B/2<W<B 


TYPE  B:  RAMP 

TYPE  D:  STREAMBANK 
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